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Abstract 

The  Theory  of  Constraints  (TOC)  is  the  foundation  for  a  computerized  scheduling 
system  called  DISASTER ™.  Although  this  system  has  proven  successful  in  many 
manufacturing  settings,  it  has  potential  limitations  due  to  the  sequential  heuristic 
process  by  which  it  schedules  constraints.  The  objective  of  this  thesis  was  to 
determine  the  extent  to  which  these  limitations  impact  the  due  date  performance  of 
schedules  created  by  DISASTER™.  This  objective  was  addressed  by  developing  an 
algorithm  to  simultaneously  schedule  multiple  constraints  in  a  job  shop  environment 
and  provide  the  optimal  schedule  for  minimized  maximum  tardiness.  This  algorithm 
was  used  to  obtain  solutions  for  a  matrix  of  job  shop  scheduling  problems,  which  were 
compared  with  solutions  obtained  by  using  DISASTER™.  This  comparison  showed 
that  DISASTER™  is  capable  of  producing  nearly  optimal  solutions  for  minimized 
maximum  tardiness,  but  that  this  capability  is  highly  dependent  on  proper  constraint 
sequencing. 
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A  COMPARISON  OF  THE  DISASTER ™  SCHEDULING  SOFTWARE 


WITH  A  SIMULTANEOUS  SCHEDULING  ALGORITHM  FOR 
MINIMIZING  MAXIMUM  TARDINESS  IN  JOB  SHOPS 


/.  Introduction 


General  Issue 

The  U  S.  Air  Force  (USAF)  has  found  the  Theory  of  Constraints  (TOC)  to  be 
useful  in  many  environments,  ranging  from  common  manufacturing  challenges  to 
battlefield  tactics.  In  the  operational  arena,  a  study  of  USAF  doctrine  dunng 
Operation  Desert  Storm  identified  mobility  as  a  constraint.  This  study  recommended 
the  use  of  TOC  to  provide  procedural  guidance  and  enhance  the  ability  of  F-16  units 
to  move  within  striking  distance  of  SCUD  missile  sights  (Cummins,  1993  :27).  The 
Air  Force  has  also  successfully  adopted  TOC  as  a  method  for  controlling  aircraft 
modification  programs  at  several  Air  Logistics  Centers  (ALCs).  For  example,  at 
Ogden  ALC  the  implementation  of  TOC  concepts  in  an  aircraft  wheel  repair  facility 
reduced  flow  days  by  75%  while  increasing  throughput  by  38%  (Demmy  and  Petnni, 
1992.6). 

In  an  ALC,  the  complex  flow  of  work  is  often  modeled  as  a  job  shop.  In  a  job 
shop,  there  may  be  any  number  of  operations  in  a  given  job,  and  the  workflow  is 
multidirectional.  As  such,  any  given  machine  may  be  required  to  process  the  job  more 
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than  once  in  its  production  sequence  (Baker,  1974:178).  Job  shop  problems  have  been 

undertaken  by  mar./  researchers  over  the  past  several  decades  with  limited  success 

The  general  job-shop  problem  is  a  fascinating  challenge  Although  it  is  easy  to 
?tate  and  to  visualize  what  is  required,  it  is  extremely  difficult  to  make  any 
progress  whatever  toward  a  solution.  Many  proficient  people  have  considered 
the  problem,  and  all  have  come  away  empty-handed.  Since  this  frustration  is 
not  reported  in  the  literature,  the  problem  continues  to  attract  investigators,  who 
just  cannot  believe  that  a  problem  so  simply  structured  can  be  so  difficult  until 
they  have  tried  it.  (Conway,  Maxwell,  and  Miller,  1967:  103) 

This  problem  continues  to  draw  attention  because  manufacturers  view 

improvements  in  production  scheduling  as  a  necessity  to  future  success  in  the  global 

market  (Wight,  1981:39)  A  1981  study  found  that  the  typical  manufacturing  job  in 

the  U  S  spent  about  ten  times  more  time  waiting  in  queue  than  in  actual  production 

Elimination  of  this  unproductive  time  could  improve  inventory  turn  times  by  a  factor 

of  one  hundred  (Kanet,  1981:58).  Clearly,  yesterday's  scheduling  methods  failed  to 

maximize  the  effectiveness  of  production  processes.  However,  since  its  introduction  in 

the  early  1 980's,  TOC  has  helped  many  job  shop  managers  to  achieve  better  system 

performance. 

TOC  simplifies  the  management  of  a  system  because  it  allows  managers  to  control 
their  processes  by  scheduling  and  monitoring  only  the  system  constraints.  In  a 
production  process,  the  application  of  TOC  usually  results  in  the  identification  of  no 
more  than  two  or  three  constraints  (Rose,  1993).  A  schedule  for  the  constraints  is 
much  simpler  than  a  schedule  for  the  entire  system,  and  may  be  obtained  using 
mathematical  solution  techniques  that  are  incapable  of  solving  more  computationally 
complex  problems. 
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The  usefulness  of  a  schedule  is  often  expressed  in  its  ability  to  minimize  or 
maximize  a  given  performance  measure  Common  performance  measures  include 
schedule  cost  and  schedule  performance  (Graves,  1981:648)  As  in  many  situations, 
cost  limitations  force  a  tradeoff  in  system  performance.  The  objective  of  a  solution 
technique  must  therefore  be  specified  (i.e.  minimize  cost,  maximize  profit,  minimize 
flow  time,  etc  )  Due-date  performance  is  often  considered  to  be  the  most  important 
schedule  performance  measure 

Of  the  various  measures  of  performance  that  have  been  considered  in  research 
on  sequencing,  certainly  the  measure  that  arouses  the  most  interest  in  those 
who  face  practical  problems  of  sequencing  is  the  satisfaction  of  preassigned  job 
due-dates.  Equipment  utilization,  work-in-process  inventory,  and  job  flow-time 
are  all  interesting  and  more  or  less  important,  but  the  ability  to  fulfill  delivery 
promises  on  time  undoubtedly  dominates  these  other  considerations.  (Conway, 
Maxwell,  and  Miller,  1967:229) 

Prominent  measures  of  due-date  performance  include:  number  of  tardy  jobs,  total 
tardiness,  mean  tardiness,  and  maximum  tardiness.  Tardiness  is  defined  as: 

tardiness  -  max  {completion  time  -  due  date,  0)  (1) 

Maximum  tardiness  has  been  chosen  as  the  topic  of  this  research,  primarily  because  of 
its  applicability  to  ALC  operations.  For  example,  if  a  number  of  jobs  must  be 
performed  in  order  to  complete  the  repair  of  an  aircraft  in  the  depot,  it  is  the  most  late 
job  that  will  determine  when  the  repairs  will  be  completed. 

As  part  of  a  strategy  to  improve  the  performance  of  their  production  and 
modification  schedules,  Wamer-Robins  ALC  and  Ogden  ALC  have  both  entered  into 
limited  partnerships  with  the  Avraham  Y.  Goldratt  Institute,  developer  of  a  TOC  based 
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computerized  scheduling  system  called  DISASTER ™.  At  Warner  Robins  ALC  alone, 
over  $500,000  has  been  spent  on  software  and  training  to  implement  DISA  STER ™ 
(Swartz,  1993). 

Since  its  introduction  in  1988,  DISASTER ™  has  allowed  many  manufacturing 
companies  to  quickly  create  schedules  that  improve  production  throughput  (Severs, 
1991:166).  However,  the  heuristic  technique  used  to  schedule  constraints  in  a 
sequential  manner  is  a  common  criticism  of  DISA  STER ™  because  a  user  may 
generate  different  schedules  for  the  same  plant  depending  on  the  order  in  which  the 
constraints  are  scheduled  (Newbold,  1992:1).  DISASTER ™  creates  a  schedule  for  one 
constraint,  and  then  builds  schedules  for  subsequent  constraints  around  the  existing 
schedule.  This  sequential  process  introduces  limitations  which  potentially  sub¬ 
optimize  the  final  schedule  (Simons,  1992:2).  The  Goldratt  Institute  acknowledges  the 
potential  for  eliminating  this  deficiency  by  simultaneously  scheduling  constraints,  but 
has  not  yet  pursued  such  an  approach  (Newbold,  1992:1).  The  USAF  also  recognizes 
the  possible  limitations  of  sequential  scheduling,  and  the  Production  Policy  Division  of 
the  Air  Force  Materiel  Command  Directorate  of  Logistics  has  sponsored  this  research 
effort  to  investigate  the  potential  for  improvement  in  this  area. 

Specific  Problem  Statement 

The  DISA  STER™  solution  algorithm  is  potentially  limited  by  the  sequential 
process  which  it  uses  to  schedule  constraints.  The  extent  to  which  these  limitations 
affect  the  maximum  tardiness  of  the  resulting  schedules  is  unknown. 
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Research  Objectives 

The  primary  goal  of  this  thesis  is  to  determine  the  extent  to  which  the  limitations 
of  sequential  scheduling  impact  the  due  date  performance  of  jobs  scheduled  by 
DISA  STER™  Our  research  objectives  are: 

1 .  Search  the  literature  to  determine  whether  an  optimal  solution  method  exists  for 
simultaneously  scheduling  multiple  constraints  in  a  job  shop. 

2.  If  an  appropriate  solution  method  cannot  be  found,  develop  a  method  to 
simultaneously  schedule  multiple  constraints  in  a  job  shop. 

3.  Compare  the  maximum  tardiness  of  schedules  created  by  DISASTER ™  with 
that  of  schedules  created  by  a  simultaneous  scheduling  method. 

4.  Determine  the  extent  to  which  the  performance  of  a  simultaneous  scheduling 
method  is  related  to  characteristics  of  the  job  shop. 

Scope 

We  will  address  our  research  problem  by  identifying  or  developing  an  algorithm  to 
simultaneously  schedule  multiple  constraints  in  accordance  with  TOC  We  will  use 
this  algorithm  to  obtain  solutions  to  job  shop  scheduling  problems,  which  we  will 
compare  with  solutions  for  the  same  problems  produced  by  DISASTER™. 

Representative  job  shop  scheduling  problems  must  be  mathematically  modeled 
before  a  solution  algorithm  can  be  applied.  A  matrix  of  representative  job  shop 
scheduling  problems  for  this  research  was  developed  by  Captain  Stewart  W.  James  and 
Captain  Bruno  A.  Mediate  in  diesis  AFIT/GS M/L AS/93  S-9.  This  matrix  consists  of 
problems  containing  two  constraints  in  a  variety  of  environments.  The  problems  were 
developed  through  the  manipulation  of  three  primary  design  factors. 
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1  Plant  type 

a.  A-plants,  where  one  finished  product  is  created  from  several  raw  materials. 

b.  T-plants,  where  several  finished  products  are  created  from  several  raw 
materials. 

c.  V-plants,  where  several  finished  products  are  created  from  one  raw 
material 

2.  Percentage  resource  criticality  factor  (%RCF)  of  the  least-loaded  constraint 
(percent  of  maximum  production  capacity  required). 

3.  The  difference  of  percentage  resource  criticality  factor  (%ARCF)  between  the 
two  constraints  (difference  in  maximum  production  capacity  required). 

Details  of  how  the  design  factors  were  combined  to  form  problems  are  included  in 

Appendix  A  and  will  be  explained  in  more  depth  in  Chapter  3. 

A  ssumptions  Several  assumptions  must  be  incorporated  in  our  algorithm  to 

account  for  the  application  of  TOC  and  to  make  our  results  comparable  to  the 

schedules  produced  by  DISA  STER™.  These  assumptions  involve  the  use  of  rods  and 

buffers,  both  of  which  will  be  explained  in  more  detail  in  the  next  chapter  The 

specific  assumptions  are: 

1.  All  buffers  are  eight  hours  long. 

2.  A  job  is  late  when  it  breaches  more  than  one-half  of  the  shipping  buffer. 

3.  All  rods  are  four  hours  long  (one  half  buffer  length). 

4.  All  constraints  are  known,  and  no  other  constraints  will  be  identified  through 
subordination. 

5.  All  operations  are  available  for  processing  any  time  after  the  ready  time 
imposed  by  the  rods  between  operations. 

In  addition  to  the  assumptions  required  to  make  a  valid  comparison  with 

DISA  STER™,  there  are  several  additional  assumptions  that  will  apply  to  our  model  of 


6 


the  job-shop.  These  assumptions  are  drawn  from  the  works  of  Conway,  Maxwell,  and 
Miller  (1967)  and  Baker  (1974). 

1.  Only  one  of  each  machine  type  exists  (i.e.,  no  parallel  processors). 

2.  All  machines  are  continuously  available  and  process  material  at  a  fixed  rate. 

3.  Each  operation  has  a  constant  processing  time  throughout  the  run  of  the  model. 

4.  Setup  times  are  either  negligible  or  included  in  processing  time  (i.e.,  setups  are 
not  sequence  dependent). 

5.  Operation  pre-emption  is  not  allowed. 

6.  Machines  are  not  interchangeable  and  can  process  only  one  job  at  a  time. 

A  pp  roach 

We  will  first  review  literature  relevant  to  TOC  and  production  scheduling, 
including  an  examination  of  the  DISA  STER ™  scheduling  program.  We  will  then 
explore  literature  pertaining  to  scheduling  problem  classification,  schedule  generation, 
and  mathematical  optimization  techniques.  In  the  methodology,  we  will  first  outline 
the  development  of  an  algorithm  to  solve  the  scheduling  problems.  We  will  then 
review  the  methods  of  data  collection  and  statistical  analysis.  The  results  of  our 
research  will  then  be  presented,  followed  by  conclusions  and  recommendations  for 
further  research. 


7 


II.  Review  of  the  Literature 


The  Theory  of  Constraints 

It  is  well  known  that  a  chain  is  only  as  strong  as  its  weakest  link.  TOC  attempts 
to  strengthen  the  chain  (the  manufacturing  process)  by  focusing  improvement  efforts 
on  the  weakest  links,  which  it  designates  as  constraints  (Simons  and  Moore,  1992:2). 

Constraints.  A  constraint  is  anything  that  restricts  the  system  from  achieving  a 
higher  level  of  performance  toward  die  organizational  goal  (Demmy  and  Petrini, 
1992:7).  Eliyahu  Goldratt,  founder  of  TOC  suggests  dealing  with  constraints  through 
a  process  of  five  focusing  steps: 

1.  Identify  the  system's  constraints). 

2.  Exploit  the  system^  constraints).  In  other  words,  determine  how  to  achieve 
the  maximum  performance  from  the  constrained  resources. 

3.  Subordinate  all  else  to  the  above  decision.  Processing  by  all  non-constraint 
resources  should  be  solely  focused  on  supporting  the  constraints). 

4.  Elevate  the  system^  constraints).  Address  the  restrictions  which  cause  die 
constraints),  and  attempt  to  eliminate  them  (i.e.  increase  capacity  of 
constrained  resource). 

5.  If  in  die  previous  steps  a  constraint  has  been  broken,  go  back  to  step  1,  bat  do 
not  allow  inertia  to  cause  a  system's  constraint  (Goldratt  and  Cox,  1992:303). 

There  are  three  basic  types  of  constraints:  external  constraints,  such  as  government 
regulations  or  market  demand;  internal  policy  constraints,  resulting  from  organizational 
procedures  and  rules;  and  internal  resource  constraints,  resulting  when  the  demand  for 
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a  resource  exceeds  its  physical  limitations  (Demmy  and  Petrini,  1992:8).  The  output 
of  one  constraint  may  later  be  the  input  for  other  constraints,  a  condition  known  as 
constraint  interaction  (Newbold,  1990.1). 

Internal  resource  constraints  can  be  classified  as  bottlenecks  or  capacity  constrained 
resources  (CCR)  (Goldratt,  1990:187).  A  bottleneck  is  a  resource  that  does  not  have 
adequate  capacity  to  satisfy  average  work  load  demands  (Goldratt,  1990:189).  A  CCR 
has  adequate  capacity  to  handle  the  average  work  load,  but  is  overwhelmed  by  work 
load  peaks  due  to  a  lack  of  protective  capacity  (Goldratt,  1990:187).  Since  it  is  only 
the  lack  of  protective  capacity  that  creates  the  capacity  constraint  for  a  CCR,  it  cannot 
be  considered  the  prime  constraint  (Goldratt,  1990:189).  Bottlenecks  will  be  the  only 
type  of  constraint  dealt  with  in  this  research  effort. 

Drum -Buffer-Rope.  The  job  shop  scheduling  technique  used  to  implement 
Goldratt's  five-step  process,  known  as  Drum-Buffer-Rope  (DBR),  is  summarized  here 
from  Goldratt  and  Fox's  1986  book.  The  Race. 

Once  a  constraint  has  been  identified,  we  must  exploit  it  and  subordinate  the  rest 
of  the  system  to  it.  This  is  where  the  drum  comes  into  action.  The  constraint  is 
scheduled  for  maximum  efficiency,  and  this  schedule  (the  drum)  dictates  the  pace  of 
the  rest  of  the  system.  Essentially,  the  constraint  is  beating  the  drum  for  the  rest  of 
the  resources  to  march  to. 

As  the  rest  of  the  resources  attempt  to  march  to  the  constraint's  drum  they  may 
encounter  problems  that  temporarily  halt  overall  production.  These  problems  are  dealt 
with  through  the  use  of  three  types  of  buffers:  constraint  buffers,  assembly  buffers. 
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and  shipping  buffers.  If  a  resource  feeding  work  to  the  constraint  is  disrupted,  the 
constraint  could  potentially  be  starved  for  work,  causing  the  entire  system  to  lose 
throughput.  A  constraint  buffer  of  material  placed  in  front  of  the  constraint  will  help 
prevent  this  occurrence.  The  size  of  the  buffer  is  measured  in  terms  of  the  amount  of 
time  it  will  take  to  process  a  given  amount  of  material.  For  example,  an  "eight-hour" 
buffer  would  consist  of  the  number  of  parts  the  constraint  could  process  in  an  eight 
hour  day.  By  providing  this  limited  supply  of  work-in-process  (WIP)  before  the 
constraint,  we  allow  it  to  continue  working  through  any  disruption  of  preceding 
resources.  Since  the  preceding  resources  have  greater  capacity  than  the  constraint  (or 
else  they  would  also  be  constraints)  they  can  use  their  excess  capacity  to  replenish  the 
buffer  after  the  disruption  has  been  corrected.  A  ssembly  buffers  are  placed  ahead  of 
assembly  processes  to  hold  non-constraint  parts  which  are  to  be  joined  with  parts  from 
a  constraint.  Assembly  buffers  ensure  that  parts  passed  from  the  constraint  continue  to 
move  through  the  production  process  without  delay  Finally,  a  shipping  buffer  is 
placed  at  the  end  of  the  assembly  process.  This  buffer  absorbs  disruptions  from  non¬ 
constraint  processes  occurring  after  processing  by  the  final  constraint  which  might 
result  in  failure  to  achieve  customer  due  dates. 

Finally,  it  is  necessary  to  control  the  release  of  raw  material  into  the  plant  at  a 
pace  that  will  not  cause  the  growth  of  excess  inventory  beyond  the  limits  of  the 
buffers.  This  is  done  by  establishing  a  control  process,  or  rope,  which  only  allows  the 
release  of  raw  materials  into  the  plant  at  a  pace  that  supports  the  constraint's  drumbeat 
(Goldratt  and  Fox,  1986:96-100). 
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DISASTER™.  The  DISASTER ™  scheduling  software  was  developed  by  the 
Avraham  Y.  Goldratt  Institute  in  1986  DISASTER™  is  a  personal-computer  based 
scheduling  software  package  that  implements  the  DBR  scheduling  technique.  The  user 
works  interactively  with  the  software  to  develop  a  schedule  for  the  constraints  in  the 
plant.  After  the  user  describes  the  layout  of  his/her  plant,  DISASTER™  determines 
which  processes  represent  constraints.  If  multiple  constraints  exist,  DISA  STER ™ 
recommends  a  sequence  in  which  to  schedule  the  constraints.  The  choice  of  which 
constraint  to  schedule  first  often  affects  the  quality  of  DISA  STER™'s  final  schedule. 
The  performance  of  DISA  STER ™  under  these  conditions  is  examined  further  in  a 
thesis  by  Captains  James  and  Mediate  (AFIT/GSM/LAS/93S-9). 

Ruins  Once  a  decision  has  been  made  on  which  constraint  to  schedule  first, 
DISASTER ™  constructs  a  diagram  of  the  ruins,  a  Gantt-chart  representation  of  the 
workload  on  the  constraints  under  the  assumption  of  infinite  capacity  (Goldratt, 
1990:220).  The  ruins  are  then  leveled  (backward  scheduled  in  reverse  order  based  on 
due  dates)  to  accommodate  the  limited  capacity  of  the  constraints  (Goldratt,  1990:204). 
If  the  quantity  of  work  exceeds  current  capacity,  this  levelling  process  will  result  in  an 
infeasible  schedule  where  work  must  begin  in  the  past  (Goldratt,  1990:205).  To 
resolve  this  problem,  the  blocks  are  then  forward  scheduled  from  the  present  to 
represent  the  best  effort  that  can  be  achieved  by  beginning  this  work  immediately 
(Goldratt,  1990:205). 

Buffers.  The  schedules  created  by  DISA  STER™  account  for  all  three  types  of 
buffers  previously  described.  DISASTER™  indicates  a  high  probability  that  a  job  will 
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be  late  if  the  final  constraint  operation  is  completed  less  than  one-half  shipping  buffer 
before  the  due  date  (Goldratt,  1990:206). 

Rods.  DISASTER™  uses  die  concept  of  rods  to  protect  the  flow  of  parts  that 
must  pass  through  multiple  constraints,  or  the  same  constraint  more  than  once 
(Goldratt,  1990:218).  Rods  prevent  two  constraint  operations  in  the  same  job  from 
being  sequenced  too  close  together.  Scheduling  constraint  operations  too  close 
together  would  increase  the  probability  that  a  system  disruption  could  starve  a 
constraint  and  delay  job  completions.  Rods  are  heuristically  defined  to  be  one-half  the 
length  of  the  constraint  buffer  and  their  placement  is  determined  by  the  processing 
time  per  part  on  the  constraints  (Goldratt,  1990:  218).  As  shown  in  Figure  1,  the 
placement  of  the  rods  can  allow  an  overlap  of  operations,  or  a  gap  between  operations 
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In  the  first  example,  two  operations  (A  and  B)  each  process  ten  parts.  The 
processing  time  per  part  for  operation  B  is  greater  than  for  operation  A.  A  forward 
rod  begins  at  the  completion  time  of  the  first  part  in  operation  A  and  extends  forward 
in  time.  A  backward  rod  begins  at  the  starting  time  of  the  last  part  in  operation  B  and 
extends  backward  in  time.  Operations  cannot  be  scheduled  any  closer  together  than  is 
permitted  by  the  rods.  In  this  example,  the  placement  of  the  rods  creates  a  gap 
between  the  completion  of  operation  A  and  the  start  of  operation  B. 

The  second  example  is  analogous,  except  that  the  predecessor  operation  (C)  has  a 
longer  processing  time  per  part  than  the  successor  operation  (D).  In  this  particular 
case,  the  placement  of  the  rods  creates  an  overlap  between  the  completion  of  operation 
C  and  the  start  of  operation  D.  It  must  be  noted  that  this  overlap  between  operations 
can  only  occur  when  the  operations  are  on  different  constraints  and  when  the  batch 
being  processed  consists  of  individual  items  which  can  be  transferred  to  the  next 
operation  independently  of  the  remaining  items  in  the  batch. 

Consistent  with  TOC  philosophy,  DISA  STER ™  does  not  provide  a  schedule  for 
non-constraints.  Non-constraints  work  on  a  first-come-first-served  basis  according  to 
the  rate  of  material  release  (a  function  of  the  rope).  DISASTER™  therefore  produces 
executable  schedules  with  significantly  less  input  (or  output)  data  than  conventional 
scheduling  packages  (Severs,  1992:167-168). 
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Schedule  Development 

Theoretically,  an  infinite  number  of  feasible  schedules  can  be  generated  for  any  job 
shop  problem,  since  an  infinitely  variable  amount  of  idle  time  can  exist  between  any 
two  adjacent  operations.  However,  when  attempting  to  maximize  any  regular  measure 
of  scheduling  performance,  it  is  clearly  desirable  to  eliminate  any  unnecessary  idle 
time  in  the  schedule  (Baker,  1974:  181). 

Semi-Active  Schedules.  Superfluous  idle  time  can  be  eliminated  by  an  adjustment 
called  a  local-left-shift:  the  movement  of  each  operation  to  the  earliest  possible  time 
without  altering  the  operation  sequences  on  any  machine  (Baker,  1974:  181)  In  this 
manner,  a  set  of  schedules  with  no  superfluous  idle  time  is  created  which  is  called  the 
set  of  semi-active  schedules  (Baker,  1974:181).  For  any  regular  performance  measure 
it  is  only  necessary  to  consider  the  set  of  semi-active  schedules,  since  the  optimal 
schedule  will  be  contained  within  this  set  (Conway,  Maxwell,  and  Miller,  1967:109) 
This  now  leaves  a  finite  set  of  schedules,  since  there  is  only  one  schedule  for  each  of 
the  finite  number  of  possible  operation  sequences  (Conway,  Maxwell,  and  Miller, 
1967:109).  Although  finite,  this  set  is  still  too  large  to  allow  complete  evaluation  for 
most  practical  problems. 

A  dive  Schedules.  A  subset  of  the  set  of  semi-active  schedules  can  be  created  by 
an  adjustment  called  a  global-left-shift  (Baker,  1974: 183).  A  global-left-shift  of  an 
operation  is  any  decrease  in  the  starting  time  of  the  operation  that  does  not  increase 
the  starting  time  of  any  other  operation  (Conway,  Maxwell,  and  Miller,  1967:111). 

This  allows  an  operation  to  shift  past  another  operation  into  an  interval  of  idle  time,  if 


14 


the  interval  is  large  enough  to  accommodate  the  shift  (Conway,  Maxwell,  and  Miller, 
1967:111).  This  subset,  known  as  the  set  of  active  schedules,  dominates  the  set  of 
semi-active  schedules;  which  means  that  it  is  sufficient  to  consider  only  the  set  of 
active  schedules  when  optimizing  any  regular  measure  of  performance  (Baker, 
1974:185). 

Schedule  Generation.  All  schedule  generation  procedures  operate  on  the  set  of 
operations  to  be  scheduled,  selecting  operations  one  at  a  time,  and  assigning  a  starting 
time  to  each  operation  (Conway,  Maxwell  and  Miller,  1967:112).  The  basic 
differences  between  schedule  generation  procedures  are  the  order  in  which  operations 
are  selected,  and  the  method  in  which  the  starting  times  are  determined  (Conway, 
Maxwell,  and  Miller,  1967:112).  These  differences  will  classify  a  schedule  generation 
procedure  as  single-pass  or  adjusting.  In  a  single-pass  schedule  generation  procedure, 
the  starting  time  of  an  operation  is  fixed  when  it  is  initially  assigned,  allowing  full 
schedule  generation  with  a  single  pass  through  the  set  of  operations  to  be  scheduled 
(Baker,  1974:187).  In  an  adjusting  schedule  generation  procedure,  the  initial 
assignment  of  an  operation  starting  time  is  tentative  and  may  be  modified  after 
subsequent  operations  have  been  scheduled  (Baker,  1974:188). 

Most  computer  programs  use  single-pass  schedule  generation  procedures  due  to  the 
difficulty  of  explicitly  stating  rules  for  adjustment  (Conway,  Maxwell,  and  Miller, 
1967:112).  This  does  not  place  any  limits  on  the  set  of  schedules  that  can  be 
generated  since  there  is  a  single-pass  procedure  capable  of  producing  any  schedule 
(Conway,  Maxwell,  and  Miller,  1967:113). 
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Schedule  generation  procedures  operate  on  the  set  of  candidate  operations:  those 
operations  which  have  no  unscheduled  predecessors  (Baker,  1974  189)  At  each  stage 
in  the  scheduling  process,  a  separate  set  of  scheduleabie  operations  exists,  along  with  a 
partial  schedule  consisting  of  all  the  operations  that  have  already  been  assigned 
starting  times  (Baker,  1974:189).  For  each  partial  schedule,  the  potential  start  time  of 
the  next  selected  operation  is  determined  by  the  completion  time  of  the  direct 
predecessor  of  that  operation,  and  the  earliest  time  the  machine  required  by  that 
operation  is  available  (Baker,  1974:189). 

A  pproaches  for  Solving  Scheduling  Problems 

Attempts  to  solve  complex  scheduling  problems  fall  into  two  general  categories: 
heuristic  solution  procedures,  and  optimal  solution  procedures. 

Heuristic  Solution  Procedures.  Heuristic  solution  methods  are  used  to  create 
feasible  schedules,  but  do  not  always  provide  the  optimal  schedule  for  a  given 
performance  measure.  Heuristics  are  commonly  applied  to  large  problems  where  the 
enormous  computational  effort  required  to  obtain  an  optimal  solution  makes  other 
approaches  impractical  (Baker,  1974: 195). 

Common  heuristic  methods  use  a  priority  rule  that  is  selected  based  on  the 
management  objective  of  the  schedule.  A  schedule  is  built  by  sequencing  the 
operations  according  to  the  selected  priority  rule.  For  the  objective  of  minimizing 
tardiness,  many  priority  rules  such  as  earliest-due-date  and  minimum-slack-time  have 
proven  useful  (Baker,  1984:195).  The  solution  technique  employed  by  DISASTER™ 
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is  heuristically  based,  and  does  not  guarantee  optimal  schedule  performance  (Newboid, 
1992:1) 

Optimal  Solution  Procedures.  Optimal  solution  procedures  mathematically  express 
the  relationships  between  the  operations  in  a  problem  and  seek  to  generate  a  schedule 
satisfying  all  of  these  relationships  while  optimizing  a  measurable  goal  (Heizer, 

Render,  and  Stair,  1993:94). 

Integer  Programming.  Integer  programming  is  a  general  purpose  approach  used 
to  optimally  solve  job  shop  scheduling  problems  (Baker,  1974:  192)  This  approach 
seeks  to  sequence  individual  operations  to  satisfy  both  sequencing  requirements  and 
noninterference  restrictions  (which  address  the  inability  of  any  machine  to  process 
more  than  one  operation  simultaneously)  (Manne,  1960:  219)  Technically,  integer 
programming  can  be  though  of  as  including  both  pure  integer  programs  and  mixed 
integer  programs.  In  pure  integer  programs,  ail  decision  variables  must  assume  an 
integer  value.  In  mixed  integer  programs,  decision  variables  are  both  continuous  and 
integer.  The  general  job  shop  scheduling  problem  is  usually  formulated  as  a  mixed 
integer  program.  A  formulation  of  the  problem  identifies  the  objective  to  be  optimized 
as  well  as  the  constraints  that  exist  to  characterize  the  problem. 

In  his  1960  paper  "On  the  Job-Shop  Scheduling  Problem,"  Manne  describes  an 
integer  programming  formulation  for  the  job  shop  problem.  This  formulation  assumes 
each  job  requires  each  machine  only  once. 

Let: 

xuk  =  the  starting  time  of  job  /  on  machine  k 
dlk  =  the  processing  time  of  job  i  on  machine  k 
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In  order  to  satisfy  the  noninterference  restrictions,  the  operations  must  precede  each 
other  by  a  sufficient  quantity  of  time  to  allow  the  completion  of  one  operation  before 
another  is  scheduled.  Accordingly,  if  there  are  two  operations  requiring  machine  k. 


then  either  , 

*  <kk 

or 


(2) 


(3) 

If  job  i  precedes  job  j,  then  xjk  must  be  dKk  days  after  xuk.  As  such,  the  integer 
programming  condition  becomes 


£  x 


(4) 


The  sequence  of  operations  is  determined  when  all  of  the  equations  in  the  problem 
formulation  are  simultaneously  solved  for  the  values  of  xul£.  For  even  small  problems 
this  set  of  constraints  can  quickly  develop  into  an  unmanageable  system  of 
inequalities.  With  only  4  machines  and  10  jobs  there  are  220  variables  and  390 
constraint  equations  (Conway,  Maxwell,  and  Miller,  1967:108).  The  increasing 
number  of  constraints  and  variables  has  forced  many  researchers  to  consider  integer 
programming  as  impractical  for  most  job  shop  problems.  (Conway,  Maxwell,  and 
Miller,  1967:109). 

Enumerative  Methods.  The  enormous  complexity  of  large  integer  programming 
problems  has  led  to  the  development  of  specialized  enumerative  approaches  to  obtain 
optimal  solutions  (Patterson,  1984:855).  These  techniques  enumerate  efficiently. 
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requiring  that  only  a  small  subset  of  the  total  number  of  possible  solutions  be 
examined  individually  (Levin  and  Kirkpatrick,  1975:349).  This  is  accomplished  by 
solving  portions  of  the  problem,  temporarily  ignoring  certain  restrictions  or  constraints. 
As  the  constraints  are  introduced  back  into  the  problem,  a  set  of  possible  schedules  is 
generated  and  ranked,  revealing  the  best  solution  (Patterson,  1984:855). 


Figure  2:  General  Branch  &  Bound  Procedure  (Simpson,  1991:14) 
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Enumerative  techniques  are  based  on  three  fundamental  concepts:  separation, 
relaxation,  and  fathoming  (Geoffnon  and  Marsten,  1972:465).  These  concepts  are 
described  below,  and  a  strategy  for  their  implementation  is  depicted  in  Figure  2. 

Separation.  Separation  (Step  1)  enables  a  "divide  and  conquer  strategy"  by 
dividing  a  problem  into  simplified  subproblems  that  are  easier  to  solve  than  the 
original  problem.  Subproblems  must  meet  the  following  two  criteria: 

1.  Every  feasible  solution  of  the  parent  problem  is  a  feasible  solution  of  exactly 
one  of  the  subproblems. 

2.  A  feasible  solution  of  any  of  the  subproblems  is  a  feasible  solution  in  the 
parent  problem  (Geoffrion  and  Marsten,  1972:466). 

The  most  common  way  of  separating  a  problem  is  by  using  contradictory  constraints 

on  a  single  variable  (Geoffrion  and  Marsten,  1972:466).  An  example  of  this  approach 

would  be  fixing  one  variable  to  a  single  value  and  solving  the  rest  of  the  problem  with 

this  variable  constant.  The  optimal  solution  for  this  subproblem  can  be  selected  (Step 

2)  as  a  candidate  for  the  optimal  solution  to  the  parent  problem  while  other 

subproblems  are  examined. 

Relaxation.  Relaxation  (Step  3)  determines  a  limit  to  the  value  of  the  objective 
function  of  a  parent  unrelaxed  problem  by  loosening  the  constraints  on  the  problem. 
This  simplification  of  the  parent  problem  implies  the  following  relationships: 

1 .  If  the  subproblem  has  no  feasible  solutions,  the  same  is  true  for  the  parent 
problem. 

2.  If  an  optimal  solution  of  the  subproblem  is  feasible  in  the  parent  problem,  then 
it  is  an  optimal  solution  of  the  parent  problem  (Geoffrion  and  Marsten, 
1972:467). 
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Relaxation  can  be  performed  by  simply  omitting  one  of  the  constraints  or  by 
eliminating  the  bounds  or  nonnegativity  conditions  on  the  variables  in  the  parent 
problem  (Geoffrion  and  Marsten,  1972:467). 

Fathoming.  Fathoming  (Step  4)  involves  discarding  subproblems  from  further 
consideration  if: 

1 .  The  subproblem  can  be  shown  to  have  no  feasible  solution. 

2.  The  subproblem  has  no  feasible  solution  better  than  the  current  candidate 
solution  to  the  parent  problem  (Geoffrion  and  Marsten,  1972:468). 

If  an  optimal  solution  to  the  subproblem  can  be  found,  then  the  subproblem  is 

completely  accounted  for,  and  there  is  no  need  for  further  separation  (Geoffrion  and 

Marsten,  1972:468). 

By  applying  the  techniques  of  separation,  relaxation  and  fathoming  to  a  problem,  it 
is  expected  that  the  solution  space  (defined  by  the  constraint  set)  can  be  systematically 
searched  to  find  and  verify  the  optimal  solution  to  the  problem  (provided  one  exists). 
One  enumerative  technique  in  which  these  principles  are  applied  is  the  branch  and 
bound  method,  where  solutions  to  valid  subproblems  provide  a  bound  to  the  optimal 
solution  of  the  initial  problem. 

Branch  and  Bound.  Branch  and  bound  (as  applied  to  a  scheduling  problem)  is  an 
enumerative  solution  method  which  creates  a  "tree"  of  partial  schedules.  Each  node  in 
the  tree  represents  an  active  schedule  at  a  stage  when  more  than  one  operation  could 
be  scheduled  on  a  particular  machine.  Each  "branch"  emanating  from  a  node 
represents  the  selection  of  one  of  the  competing  operations  (Conway,  Maxwell,  and 


21 


Miller,  1967:117).  Each  branching  operation  can  create  as  many  new  partial  schedules 
as  there  are  competing  operations  at  that  stage.  Creation  of  all  possible  branches  at 
each  node  will  eventually  result  in  a  set  of  all  possible  feasible  schedules  (Stinson, 
Davis,  and  Khumawala,  1978:253).  This  set  is  too  large  to  be  practical  for  finding  the 
optimal  solution  to  even  medium-sized  problems.  Therefore,  partial  schedules  are 
fathomed  using  lower  bounds  and  upper  bounds  as  a  means  to  reduce  the  time 
required  to  find  and  verify  the  optimal  solution. 

A  lower  bound  based  upon  an  existing  partial  schedule  is  calculated  for  every 
node,  and  represents  the  "best  case"  value  for  the  objective  function  for  all  of  the 
complete  schedules  which  ultimately  emanate  from  that  node  (Conway,  Maxwell,  and 
Miller,  1967:117).  The  value  of  a  lower  bound  may  not  be  feasible  in  a  complete 
schedule,  because  it  is  often  calculated  through  relaxation  of  the  constraints.  An  upper 
bound  is  calculated  after  a  complete  schedule  has  been  generated.  It  consists  of  the 
value  of  the  objective  function  measure  for  the  complete  schedule  (Conway,  Maxwell 
and  Miller,  1967:  117). 

Techniques  known  as  search  methods  use  these  bounds  to  converge  on  an  optimal 
solution  by  evaluating  the  bounds  of  the  branches  and  determining  how  the  branch- 
and-bound  tree  is  developed  (Johnson,  1988:17).  An  example  of  a  search  method  is 
the  last-in-first-out  search  (LIFO).  This  method  begins  with  the  initial  node  in  the 
solution  tree,  forms  all  branches  descending  from  that  node,  and  selects  the  branch 
with  the  lowest  lower  bound.  A  node  is  then  formed  from  this  branch,  and  the 
branching  process  is  repeated.  An  upper  bound  (incumbent  solution's  objective 
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function  value)  is  calculated  when  all  operations  have  been  sequenced.  The  LIFO 
search  then  backtracks  along  the  node  path  to  the  nearest  node  that  has  a  descending 
node  that  has  not  been  extended.  If  this  descending  node's  lower  bound  is  superior  to 
the  upper  bound,  the  branching  process  resumes  (Johnson,  1988:24).  Branches  with 
lower  bounds  inferior  to  die  upper  bound  are  fathomed  (i.e.  eliminated  from  further 
consideration).  Each  time  a  new  complete  schedule  is  produced,  its  objective  function 
value  is  compared  to  the  current  upper  bound.  If  this  new  value  is  superior,  the  upper 
bound  is  adjusted  and  the  schedule  which  produced  it  becomes  the  incumbent.  If  an 
optimal  solution  exists,  this  process  will  continue  until  an  incumbent  solution 
backtracks  all  the  way  to  the  initial  node,  and  is  thus  identified  as  the  optimal  solution 

The  value  of  fathoming  lies  in  the  fact  that  it  reduces  the  number  of  alternative 
schedules  which  must  be  explicitly  considered.  A  node  whose  lower  bound  exceeds 
the  current  upper  bound  cannot  possibly  produce  an  optimal  solution.  We  know  this 
because  the  existence  of  the  upper  bound  indicates  that  we  already  have  a  schedule 
better  than  the  best  possible  case  which  could  be  obtained  from  the  node  being 
fathomed. 

Graph-Theoretical  Method.  One  method  of  implementing  die  branch-and-bound 
solution  technique  for  the  job  shop  scheduling  problem  is  the  Graph-Theoretical 
machine  sequencing  algorithm  outlined  by  Florian,  Trepant,  and  McMahon  in  their 
1971  paper  "An  Implicit  Enumeration  Algorithm  for  the  Machine  Sequencing 
Problem."  Their  technique  involves  a  single-pass  method  which  generates  all  active 
schedules,  using  a  LIFO  search  method  to  obtain  an  optimal  solution  for  regular 
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performance  measures.  The  following  sections  summarize  the  portions  of  their 
formulation  and  solution  technique  relevant  to  this  research  effort. 

A  graphical  notation  is  used  to  represent  precedence  conditions  between  operations. 
This  notation  creates  a  node  for  every  operation,  along  with  an  initial  node  (0),  and  a 
terminal  node  (*).  Conjunctive  arcs  are  one-way  arrows  which  connect  nodes 
representing  successive  operations  within  each  job.  These  arcs  reflect  precedence 
constraints  between  operations.  Machine  interference  constraints  that  must  be  resolved 
to  create  a  complete  schedule  are  represented  by  opposing  pairs  of  disjunctive  arcs. 
Table  1  provides  data  for  a  sample  problem  consisting  of  two  machines  and  three  jobs. 


Table  1.  Sample  Problem  Data  (Florian,  Trepant,  and  McMahon,  1971  :B783) 


Operation 

Job 

[2222 

Duration  j 

1 

1 

i 

4 

2 

1 

2 

3 

3 

2 

2 

2 

4 

2 

1 

3 

5 

3 

1 

1 

6 

3 

2 

2 

Figure  3  depicts  this  problem  in  graph-theoretical  form.  Note  that  conjunctive  arcs 
are  shown  between  successive  operations  in  the  same  job  (such  as  operations  1  and  2). 
By  contrast,  disjunctive  arcs  join  operations  which  must  be  processed  on  the  same 
machine  (such  as  operations  1  and  4),  representing  machine  interference  constraints. 

A  solution  for  the  problem  is  obtained  by  determining  the  optimal  sequence  of 
operations  on  each  machine.  This  is  sequence  is  created  by  "resolving"  the  disjunctive 
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Figure  3:  Graph-Theoretical  Network  (Florian,  Trepant,  and  McMahon, 


1971  B783) 

arcs  and  replacing  them  with  conjunctive  arcs,  indicating  which  of  the  two  operations 
will  be  scheduled  earlier. 

Resolving  the  Disjunctive  A  res.  A  disjunctive  arc  is  resolved  when  the 
sequence  of  the  two  operations  connected  by  the  arc  is  determined.  An  operation  is 
scheduled  on  a  machine  when  all  disjunctive  arcs  connected  to  the  operation  have  been 
resolved.  To  begin  this  process,  it  is  necessary  to  identify  the  initial  cut  (CJ.  A  cut 
is  defined  as  the  set  of  all  candidate  operations  (those  with  no  unscheduled 
predecessors).  In  this  example,  C0  consists  of  operations  1,  3,  and  5,  the  only 
operations  with  no  predecessors  other  than  die  initial  node.  A  generating  set  is 
obtained  from  this  cut,  consisting  of  all  operations  in  die  cut  on  the  same  machine  as 
the  operation  in  the  cut  with  the  earliest  possible  finish  time.  A  branch  and  bound  tree 
node  is  created  for  each  alternative  within  the  generating  set.  Examination  of  Table  1 
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reveals  that  operation  5  has  a  duration  of  1  time  unit.  Since  all  operations  in  the 
initial  cut  are  able  to  start  at  time  zero,  operation  5  has  the  earliest  possible  finish 
time.  The  generating  set  consists  of  all  operations  in  C0  that  are  on  the  same  machine 
as  operation  5  (  machine  1 ).  Thus  we  see  that  the  generating  set  for  C0  is  {7.5}. 

Since  there  is  more  than  one  operation  in  the  generating  set,  it  is  necessary  to 
decide  which  to  sequence  first  on  machine  1 .  A  decision  is  made  to  schedule  one  of 
the  nodes  in  the  generating  cut  by  calculating  and  comparing  the  lower  bounds  of  the 
competing  operations.  The  operation  with  the  lowest  bound  is  selected  and  scheduled 
by  resolving  all  disjunctive  arcs  connected  to  the  node.  Because  the  lower-bound 
calculation  method  used  by  Florian,  Trepant,  and  McMahon  is  not  relevant  to  this 
research  effort,  it  will  not  be  reviewed  here,  and  we  will  simply  assume  that  operation 
7  has  the  lowest  bound. 

Since  operation  7  has  the  lowest  bound,  it  becomes  the  first  operation  scheduled  on 
machine  1 .  A  new  cut  is  now  defined  to  begin  selection  of  the  next  operation  to 
schedule.  This  new  cut  (C,)  is  created  by  removing  operation  7  from  the  cut  C0,  and 
adding  its  successor,  operation  2  to  the  cut.  A  generating  set  is  obtained  from  this  cut 
as  before,  and  the  procedure  repeats  until  all  operations  have  been  scheduled.  This 
method  is  known  as  the  consecutive  cut  enumeration  method.  As  this  method 
progresses,  a  tree  of  nodes  begins  to  develop,  such  as  in  Figure  4.  Florian,  Trepant, 
and  McMahon  have  proven  that  this  method  generates  all  active  schedules  for  the  job 
shop  problem  (Florian,  Trepant,  and  McMahon,  1971  :B-786).  As  mentioned 
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previously,  it  is  sufficient  to  consider  only  the  set  of  active  schedules  when  seeking 
the  optimal  solution  for  any  regular  performance  measure  (Baker,  1974:185) 


Figure  4:  Tree  of  Nodes 


A  complete  schedule  is  generated  when  all  disjunctive  arcs  have  been  resolved. 
Figure  5  contains  a  network  representation  of  a  completed  schedule  that  maintains  the 
precedence  and  machine  relationships  for  all  operations  in  our  example  problem. 


Figure  5:  Resolved  Graph-Theoretical  Network 


A  Gantt  chart  corresponding  to  the  schedule  shown  in  Figure  S  is  provided  in  Figure  6 
to  aid  in  the  interpretation  of  the  network. 


Figure  6:  Gantt  Chart  for  Resolved  Network 


Backtracking.  Once  a  complete  schedule  has  been  generated,  a  value  for  the 
objective  function  is  computed.  (Florian,  Trepant,  and  McMahon's  objective  function  is 
to  minimize  makespan).  This  becomes  the  incumbent  solution  (upper  bound).  The 
optimality  of  this  solution  is  then  tested  by  backtracking  through  the  tree  of  nodes 
The  lower  bound  of  all  nodes  in  the  generating  set  of  each  cut  are  examined.  If  any 
of  the  lower  bounds  are  less  than  the  incumbent  upper  bound,  then  the  corresponding 
operation  is  substituted  for  the  operation  originally  scheduled  from  that  cut.  A  new 
series  of  cuts  is  then  defined,  and  the  procedure  repeats  until  we  have  backtracked  all 
the  way  to  the  initial  cut.  If,  however,  a  node's  lower  bound  is  greater  than  the 
incumbent  upper  bound,  then  the  node  is  fathomed  and  the  schedules  which  would 
have  been  produced  from  that  node  need  never  be  explicitly  considered.  The  sequence 
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producing  the  final  upper  bound  is  an  optimal  sequence  for  the  desired  management 
objective  (Florian,  Trepant,  and  McMahon,  1971  :B-782  -  B-790). 

Summary 

Our  overview  of  TOC  has  reviewed  the  rationale  allowing  die  scheduler  to  focus 
on  the  constraints  in  the  system.  Our  review  of  the  DISASTER ™  software  package 
demonstrated  how  this  program  applies  TOC  to  generate  a  useable  schedule.  We  also 
noted  the  potential  shortcoming  of  DISASTER ™  resulting  from  its  heuristic  solution 
method.  We  reviewed  the  classification  of  production  scheduling  problems,  and  the 
methods  that  exist  to  generate  schedules  for  these  problems.  Finally,  we  examined 
methods  for  determining  the  solution  to  a  production  scheduling  problem,  with 
particular  emphasis  on  the  characteristics  of  the  branch-and-bound  solution  technique 
and  the  Graph-Theoretical  method  for  producing  optimal  solutions. 
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Ill  Methodology 


In  our  survey  of  the  existing  literature,  we  did  not  find  an  optimal  solution  method 
that  met  our  requirements.  In  accordance  with  our  second  research  objective,  we 
developed  a  solution  method  for  simultaneously  scheduling  multiple  constraints  m  a 
job  shop.  The  first  step  in  developing  this  method  was  to  formulate  the  problem.  Our 
formulation  seeks  to  optimize  due  date  performance  by  minimizing  the  maximum 
tardiness  of  all  jobs.  The  algorithm  we  developed  to  solve  this  problem  was  an 
adaptation  of  the  Graph-Theoretical  method  developed  by  Floiian,  Trepant,  and 
McMahon  for  the  general  job  shop.  The  basic  node  expansion  logic  of  their  method 
was  retained,  but  new  bounds  were  required  due  to  the  different  objective  function 
Additional  modifications  were  required  to  accommodate  the  differences  between  the 
general  job  shop  and  a  job  shop  applying  the  principles  of  TOC. 

Problem  Formulation 

Our  problem  formulation  was  adapted  from  Baker's  1974  integer  programming 
formulation  for  the  job  shop  problem.  Our  formulation  uses  the  following  notation. 
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Table  2:  Function  Notation 


n 

=  the  number  of  jobs 

<7, 

=  the  number  of  operations  in  job  i 

m 

=  the  number  of  machines 

i 

=  l../»jobs 

j 

=  l..q,  operations 

k 

=  1  ..m  machines 

u 

=  completion  time  of  operation  j  of  job  /  on  machine  k 

A, 

=  due  date  of  job  / 

=  processing  time  of  operation  q  of  job  p  on  machine  k 

8tjk 

=  the  ready  time  of  operation  j  of  job  i  on  machine  k 

X 

=  very  large  positive  number 

c 

=  the  set  of  all  jobs 

=  1  if  operation  j  of  job  /  precedes  operation  q  of  job  p  on  machine  k, 
else  =  0 

(5) 

(6) 

(7) 

(8) 
(9) 


The  objective  is  to  minimize  the  maximum  tardiness  (the  difference  between  due 
date  and  completion)  across  the  set  of  all  jobs.  Equation  (7)  states  that  no  operation 


can  start  any  earlier  than  its  ready  time.  Equations  (8)  and  (9)  prevent  the  scheduling 
of  two  operations  on  the  same  machine  simultaneously.  Equation  (10)  prevents 
scheduling  of  operations  in  the  past.  The  gl]k  are  used  to  accommodate  the  rods. 


Algorithm  Development 

Our  algorithm  was  based  on  the  Graph-Theoretical  machine  sequencing  algorithm 
developed  by  Florian,  Trepant,  and  McMahon  described  in  the  previous  chapter.  This 
method  was  selected  because  their  solution  technique  was  easily  adaptable  to  the 
objective  of  minimizing  maximum  tardiness  by  simply  using  a  different  lower  bound. 
The  bound  we  developed  for  this  research  effort  is  calculated  by  assigning  each 
operation  the  due  date  of  its  job  and  sequencing  the  operations  by  earliest  due  date 
(EDD)  with  ties  broken  arbitrarily.  The  value  of  this  lower  bound  reflects  the 
minimum  possible  maximum  tardiness  at  each  point  in  sequence  of  operations.  A 
proof  of  the  validity  of  this  lower  bound  is  contained  in  the  following  sections. 

A  ssumptions.  All  operations  are  ready  for  processing  at  t=0  (referred  to  later  as  the 
assumption  of  readiness). 


Notation : 
C 
ij 
B, 
f 
h, 

T, 

4 


=  the  set  of  all  operations 
=  two  operations  within  C 

=  the  set  of  operations  preceding  operation  i  in  the  current  schedule 
=  the  completion  time  of  operation  i 
=  the  due  date  of  the  job  containing  operation  i 
=  the  tardiness  of  operation  i 
=  the  processing  time  of  operation  / 
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Single  Machine  Tardiness  Determination  Hypothesis :  On  a  single  machine,  for  all 
jobs,  tardiness  of  a  job  will  be  determined  by  the  maximum  tardiness  of  the  operations 
in  that  job. 

Single  Machine  Tardiness  Determination  Proof.  Each  job  consists  of  multiple 
operations  and  is  assigned  a  job  due  date.  Assume  that  there  exists  a  job,  A , 
consisting  of  two  operations,  i  and  j.  Job  A  will  be  complete  upon  the  completion  of 
operations  /  and  j.  It  can  therefore  be  stated  that  the  due  date  for  the  completion  of 
operations  /  and  j  is  the  due  date  for  the  completion  of  job  A .  Tardiness  is  defined  as 
the  positive  difference  between  the  completion  of  an  operation  or  job  and  its  due  date 
Since  the  due  date  is  the  same  for  both  operations  z  and  j,  the  tardiness  of  each 
operation  will  be  determined  by  the  operation's  finish  time.  In  the  two-operation 
example  with  job  A ,  there  are  only  four  possible  outcomes: 

1 .  If  /,  <  h,  and  f  <  hi  then  TA  =  0 

2.  Iff  >f  and/,  >/z„  then  TA  =  T, 

3.  If  /  >/  and  /  >hp  then  TA  =  f 

4  If/  =/  and /  >  ht  then  TA  =  Tt  =  T} 

The  above  example  can  be  extended  to  jobs  consisting  of  more  than  two  operations; 

however,  in  all  cases  the  tardiness  of  a  job  is  determined  by  the  tardiness  of  the  most 

tardy  operation. 

Single  Machine  Tardiness  Minimization  Hypothesis :  Since  the  tardiness  of  a  job 
will  be  determined  by  the  tardiness  of  its  most  tardy  operation,  the  minimized 
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maximum  tardiness  of  multiple  jobs  can  be  determined  by  minimizing  maximum 
tardiness  of  the  operations  of  these  jobs. 

Single  Machine  Tardiness  Minimization  Proof  .  The  following  proof  of  minimizing 
maximum  tardiness  with  EDD  sequencing  is  adapted  from  Baker  (1974). 

An  initial  sequence,  S,  may  be  formed  by  scheduling  the  operations  in  any  order. 
The  assumption  of  readiness  results  in  a  schedule  with  no  idle  time  •  a  continuous 
period  of  activity  on  the  machine  from  t  =  0  until  the  completion  of  the  last  scheduled 
operation.  The  makespan  of  this  activity  on  the  machine  will  remain  the  same 
regardless  of  the  sequence  of  operations. 

If  the  operations  are  not  in  EDD  sequence,  there  will  exist  two  operations,  ;  and  j, 
with  j  directly  following  i  such  that  /»,  >  hy  A  new  sequence.  S',  can  be  constructed  in 
which  /  and  j  are  interchanged  and  all  other  operations  complete  in  the  same  time  as  in 
sequence  S. 

The  tardiness  of  the  two  operations  in  sequence  S  will  be  determined  as  the 
positive  difference  between  the  completion  of  an  operation  and  its  due  date: 

T£S)  =  maxifm+drhit  0)  (10) 


Tj(S)  =  maxtfn+dsdj-hp  0) 


(ID 


The  tardiness  of  the  two  operations  in  sequence  S'  will  be  similarly  determined  as: 

Ti(S/ )  =  maxifa+dj+di-hp  0)  (12) 

TjLS')  =  max{fm+dj-hj,  0)  (13) 

If  neither  deration  has  tardiness  >  0,  there  is  no  benefit  or  penalty  in  terms  of 
minimized  maximum  tardiness  by  interchanging  the  two  operations.  If  both  operations 
have  positive  tardiness,  the  above  relationships  show  that: 

TjiS)>T£S')  (14) 

and 

r/s)  >  r/s' )  (is) 

The  tardiness  of  operation  j  in  sequence  S  is  greater  than  the  tardiness  of  /  or  j  in 
sequence  S'.  If  operation  i  or  operation  j  are  not  the  most  tardy  operations  in  the 
sequence  S,  the  value  for  Tmm  will  be  unchanged  in  sequence  S'.  This  procedure  may 
be  repeated  for  all  operations  in  the  sequence  where  the  succeeding  operation's  due 
date  is  earlier  than  that  of  the  preceding  operation.  At  each  iteration,  the  value  of  Tmm 
will  remain  the  same  or  be  reduced.  After  all  of  these  cases  have  been  resolved,  the 
resulting  sequence  will  minimize  maximum  tardiness  for  all  operations  in  the  schedule. 

Multiple  Machine  Hypothesis :  The  procedure  that  we  have  just  reviewed  for 
sequencing  operations  on  a  single  machine  can  be  used  as  a  bound  in  the  multiple 
machine  case  by  ignoring  operational  precedence  constraints  and  scheduling  each 
machine  independently.  By  scheduling  the  remaining  operations  of  all  jobs  on  each 
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applicable  machine  by  EDD,  we  will  minimize  the  maximum  tardiness  of  each  job 
relative  to  each  machine.  The  value  for  the  bound  will  then  be  the  maximum  tardiness 
of  all  jobs  among  all  machines. 

Multiple  Machine  Proof.  As  presented  in  Geoffnon  and  Marsten  (1972),  any 
constrained  optimization  problem,  P,  can  be  relaxed  by  loosening  its  constraints, 
resulting  in  a  new  problem,  P'.  In  the  case  of  our  bound,  the  problem  of  sequencing 
the  operations  among  both  machines,  P,  is  simplified  by  relaxing  the  operational 
precedence  constraints  and  the  machine  interference  constraints.  The  new  problem,  P\ 
requires  that  we  minimize  the  maximum  tardiness  among  each  machine  independently 
when  all  operations  are  available  for  processing  at  T=0.  The  property  of  relaxation 
implies  that  the  minimal  value  of  our  performance  measure  for  P  is  no  less  than  the 
minimal  value  of  the  performance  measure  for  P'.  Additionally,  if  an  optimal  solution 
of  P'  is  feasible  in  P,  then  it  is  an  optimal  solution  for  P. 

Statement  of  the  Basic  A  Igorithm 

Our  algorithm,  as  adapted  from  die  algorithm  oudined  in  Florian,  Trepant,  and 
McMahon's  1971  paper  to  reflect  the  objective  of  minimizing  maximum  tardiness,  is 
presented  below. 

1.  Problem  initialization 

a.  Designate  the  initial  cut,  a  set  consisting  of  all  operations  whose  only 
predecessor  is  the  initial  node. 

b.  Compute  earliest  start  times  for  all  operations. 
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c.  Set  upper  bound  =  oo. 

2.  Extract  Generating  Set  from  Cut 

a.  Determine  the  operation  in  the  cut  with  the  earliest  finish  time. 

b.  All  operations  in  the  cut  on  the  same  machine  as  the  above  operation  are  in 
the  generating  set. 

c.  If  there  are  no  operations  in  the  generating  set,  then  go  to  step  8. 

d.  If  there  is  only  one  operation  in  the  generating  set,  then  go  to  step  6 

e.  If  there  is  more  than  one  operation  in  the  generating  set,  then  continue. 

3.  Compute  the  lower  bound  for  every  operation  in  the  generating  set. 

a.  Assume  the  selected  operation  is  scheduled. 

b.  Sequence  all  unscheduled  operations  by  earliest  due  date. 

c.  Calculate  the  tardiness  for  each  job. 

d.  Record  maximum  tardiness.  This  is  the  lower  bound  for  the  selected 
operation. 

4.  Select  the  operation  in  the  generating  set  with  the  lowest  value  for  the  lower 
bound.  If  this  lower  bound  >  the  incumbent  upper  bound,  then  go  to  step  7; 
otherwise  schedule  the  selected  operation  next  on  its  designated  machine. 

5.  Set  the  lower  bound  for  the  selected  operation  =  oo. 

6.  Remove  the  selected  operation  from  the  cut,  and  add  any  successors  to  die 
selected  operation  to  the  cut.  Compute  new  earliest  start  and  finish  times  for 
each  unscheduled  operation.  Go  to  step  2. 

7.  If  the  current  cut  is  the  initial  cut,  then  go  to  step  9;  otherwise  backtrack  to  the 
previous  cut  and  go  to  step  4. 

8.  A  complete  sequence  has  been  generated.  Compare  the  maximum  tardiness  of 
this  sequence  to  the  incumbent  upper  bound.  If  this  sequence  results  in  a  lower 
maximum  tardiness,  then  the  maximum  tardiness  resulting  from  this  sequence  is 
the  new  upper  bound.  Go  to  step  7. 
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9.  The  last  complete  sequence  generated  was  the  optimal  sequence  for  minimizing 
maximum  tardiness. 


figure  7:  Algorithm  Flowchart  (adapted  from  Florian,  Trepant,  and  McMahon, 


1971:B-788) 
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A  daptation  of  the  Basic  A  Igorithm  to  a  TOC  Job  Shop 

Three  basic  changes  were  made  to  adapt  the  previously  defined  algorithm  to  a 
TOC  job  shop  and  allow  a  valid  comparison  to  results  generated  by  DISA  STER ™ 
Constraints.  DISASTER™  only  creates  a  schedule  for  constraint  resources,  and  it 
controls  all  non-constraint  resources  through  subordination.  In  order  to  provide 
schedules  comparable  with  those  generated  by  DISASTER ™  we  ignored  all  non¬ 
constraint  operations,  and  built  a  schedule  for  the  constraint  operations  only. 
Subordination  of  non-constraint  operations  to  our  constraint  schedule  is  beyond  the 
scope  of  this  research  effort. 

Shipping  Buffer.  DISASTER™  defines  a  job  as  tardy  if  the  last  constraint 
operation  is  completed  less  than  one-half  of  a  shipping  buffer  before  the  due  date  of 
the  job.  If  the  operation  is  completed  after  this  point,  DISASTER™  adds  one  full 
shipping  buffer  to  the  completion  time  of  the  last  constraint  operation  to  arrive  at  the 
anticipated  completion  time  of  the  job.  To  allow  a  comparison  with  results  generated 
by  DISASTER™,  we  have  used  the  same  method  to  compute  completion  times. 

Rods.  Since  we  are  only  scheduling  constraints,  TOC  mandates  the  use  of  rods 
between  operations.  In  order  to  introduce  this  concept  into  our  algorithm,  we  calculate 
an  adjusted  ready  time  for  each  operation.  This  new  ready  time  is  the  earliest  time 
that  an  operation  can  be  scheduled  without  violating  a  rod.  To  allow  a  comparison 
with  results  generated  by  DISASTER™,  we  used  the  same  half-buffer  rod  length,  and 
the  same  method  of  calculating  rod  placement. 
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Model  Coding 

The  algorithm  was  encoded  using  Borland  Turbo  Pascal  for  Windows  1.5.  The 
program  was  compiled  and  run  on  a  Unisys  personal  computer  (IBM  compatible)  with 
a  386DX  3  3 MHz  processor.  This  machine  had  eight  mega-bytes  of  RAM  and  a  387 
coprocessor.  The  operating  environment  was  Microsoft  DOS  5.0  with  Microsoft 
Windows  3.1.  A  copy  of  the  program  file  is  included  in  Appendix  B. 

Experimental  Factors 

In  accordance  with  our  fourth  research  objective,  a  split-plot  experimental  design 
was  created  to  manipulate  several  job  shop  characteristics.  The  performance  of  our 
algorithm  was  examined  through  the  manipulation  of  three  factors  to  observe  the  effect 
on  a  set  of  dependent  variables.  This  experimental  design  was  developed  by  Captain 
Stewart  W.  James  and  Captain  Bruno  A.  Mediate  in  thesis  AFIT/GSM/LAS/93S-9  to 
produce  a  set  of  benchmark  problems  and  DISASTER ™  solutions.  A  description  of 
each  factor  and  its  levels  are  provided  below. 

Plait  Type.  The  first  experimental  factor  differentiates  among  three  basic  plant 
types  used  in  this  research  effort:  A -plants,  T -plants,  and  V -plants.  A  single 
configuration  of  each  plant  type  was  used  in  the  benchmark  problems. 

An  A  -plait  is  a  plant  in  which  operations  converge,  with  a  few  finished  products 
created  from  many  raw  materials  (Fawcett  and  Pearson,  1991:50).  For  the  purpose  of 
this  research  effort,  an  A-plant  will  create  one  finished  good  from  nine  raw  materials. 

A  diagram  of  the  A-plant  used  in  this  research  is  provided  in  Appendix  A. 
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A  T-plant  is  a  plant  in  which  several  finished  goods  are  created  from  several  raw 
materials.  This  type  of  plant  is  usually  an  assemble-to-order  operation  with  a  number 
of  common  parts  used  to  assemble  the  finished  products  (Fawcett  and  Pearson, 
1991:51).  For  the  purpose  of  this  research  effort,  a  T-plant  will  create  five  finished 
goods  from  four  raw  materials.  A  diagram  of  the  T-plant  used  in  this  research  is 
provided  in  Appendix  A 

A  V -plant  is  a  plant  in  which  operations  diverge,  with  many  finished  products 
created  from  a  few  raw  materials  (Fawcett  and  Pearson,  1991:50).  For  the  purpose  of 
this  research  effort,  a  V-plant  will  create  five  finished  goods  from  one  raw  material. 

A  diagram  of  the  V-plant  used  in  this  research  is  provided  in  Appendix  A. 

Plants  exhibiting  characteristics  of  more  than  one  of  these  plant  types  are  called 
combination  plants  (Fawcett  and  Pearson,  1991:50).  Combination  plants  will  not  be 
addressed  in  this  research  effort. 

Percent  Resource  Criticality  Factor.  Of  the  ten  resource  types  used  in  the 
previously  defined  plants.  Captains  James  and  Mediate  defined  two  constraints,  the 
Blue  resource  and  the  Gold  resource.  They  operationalized  the  demand  for  each 
constraint  as  a  percent  resource  criticality  factor  (%RCF).  %RCF  is  defined  as  the 
ratio  of  total  time  a  resource  is  required  for  a  given  time  period  to  the  total  capacity 
available  for  that  resource  (Gargeya,  1992:  3).  The  second  experimental  factor 
expresses  demand  for  the  lower-constrained  resource  in  terms  of  %RCF  at  three 
different  levels:  105%,  115%,  and  125%.  These  levels  of  %RCF  are  achieved  in  the 
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benchmark  problems  by  varying  the  required  processing  time  of  each  operation  on  the 
constraints,  as  depicted  in  Appendix  A. 

Percent  Delta  Resource  Criticality  Factor.  Differences  in  demand  between  the  two 
constraint  resources  is  expressed  in  terms  of  their  third  experimental  factor  (%ARCF) 
This  factor  is  expressed  at  three  different  levels:  0%,  25%,  and  50%.  These  levels  of 
%ARCF  are  achieved  in  the  benchmark  problems  by  varying  the  relative  processing 
time  of  each  operation  on  the  constraints,  as  depicted  in  Appendix  A. 

Background  Variables.  Captains  James  and  Mediate  identified  sixteen  background 
variables,  all  but  two  of  which  were  controlled  by  holding  diem  constant.  The 
exceptions  were  product  types,  and  job  order  due  dates. 

Product  Types.  To  achieve  a  wide  range  of  conditions,  a  series  of  product 
types  were  defined,  using  various  constraint  process  sequences.  The  following  table 
identifies  the  constraint  sequence  for  each  plant  type's  products. 


Table  3:  Product  Type  Constraint  Sequence  (James  and  Mediate,  1993:115-141) 


Plant  Type 

Product  Type 

Constraint  Sequence 

A 

FG-D  Jl 

Blue-Gold-Blue 

T 

FG-B 

Blue-Blue-Gold-Blue 

T 

FG-C 

Blue-Blue-Gold 

T 

FG-D 

Gold 

T 

FG-F 

Gold 

T 

FG-G 

Blue-Blue-Gold-Gold  | 

V 

FG-A 

Blue-Blue  | 

V 

FG-B 

Blue-Gold 

V 

FG-D 

Blue 

V 

FG-E 

Gold-Gold 

V 

FG-G 

Gold-Gold-Blue 
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Due  Dates.  Job  order  due  dates  were  assigned  according  to  a  due  date 
assignment  rule  which  is  based  on  the  job  order's  arrival  date,  and  the  number  of 
processing  stations  for  the  corresponding  product  type  (Ragatz  and  Mabert,  1984:29) 
There  were  ten  job  orders  for  each  replication.  The  product  types  were  drawn  from 
the  table  above,  and  the  arrival  dates  were  randomly  selected  without  replacement 
from  a  two  week  period.  This  process  was  used  to  generate  four  replications  for  each 
treatment  of  the  experimental  factors,  resulting  in  108  benchmark  problems  (3  plant 
types  x  3  %RCF  levels  x  3  %ARCF  levels  x  4  replications  =  108  replications).  The 
replications  for  each  plant  type  were  not  identical,  hence  the  split-plot  experimental 
design.  Details  on  the  levels  of  the  experimental  factors  for  each  replication  are 
included  in  Appendix  A.  Since  the  problems  involved  two  constraints.  Captains  James 
and  Mediate  ran  each  replication  through  DISASTER™  twice,  reversing  the  sequence 
of  the  constraint  schedule  and  thereby  producing  two  separate  schedules  for  each 
replication.  Because  of  the  simultaneous  scheduling  nature  of  our  algorithm,  we 
produced  only  one  schedule  for  each  replication. 

Data  Collection 

We  used  the  108  problems  developed  by  Captains  James  and  Mediate  to  create  an 
input  data  file  for  our  algorithm.  Running  our  program  with  this  input  file,  we 
attempted  to  produce  an  optimal  solution  for  each  problem.  We  discovered  that 
several  problems  would  not  produce  a  solution  within  a  reasonable  amount  of 
processing  time.  As  a  result,  we  terminated  processing  for  each  problem  after  S 
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minutes  of  CPU  time  if  an  optimal  solution  had  not  yet  been  verified.  At  that  point, 
we  designated  the  incumbent  solution  as  a  "heuristic"  solution.  These  solutions  may, 
in  fact,  have  been  optimal,  but  the  time  required  to  verify  the  optimality  of  the 
solution  exceeded  the  designated  time  limit. 

Statistical  A  nalysis 

In  order  to  satisfy  research  objectives  three  and  four,  we  performed  a  series  of 
statistical  comparisons  using  descriptive  statistics  and  analysis  of  variance  (ANOVA) 
All  analyses  were  conducted  using  the  STAT1STIX ™  4.0  software  package 

Statistical  Comparisons.  In  order  to  compare  the  maximum  tardiness  of  schedules 
created  by  DISA  STER ™  with  that  of  schedules  created  by  a  simultaneous  scheduling 
method  (research  objective  three),  we  must  compare  the  samples  of  data  collected 
using  both  solution  methods.  Tests  for  statistical  comparisons  among  the  means  or 
medians  of  two  samples  of  data  (i.e.  paired  t-test,  sign  test)  require  independence  of 
the  data  within  each  sample  set.  Our  split-plot  experimental  design  introduces 
dependence  among  these  values,  since  the  same  four  replications  are  repeated  within 
each  setting  for  plant  type  (this  is  a  characteristic  of  the  split  plot  design)  Therefore, 
statistical  tests  of  the  means  or  medians  of  these  samples  are  not  valid  (Reynolds, 
1993).  As  a  result,  we  analyzed  all  comparisons  between  die  two  solution  methods  of 
data  qualitatively. 

Analysis  of  Variance.  Our  fourth  research  objective  requires  an  investigation  of 
the  extent  to  which  the  performance  of  a  simultaneous  scheduling  method  is  related  to 
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the  characteristics  of  the  job  shop  (research  objective  four),  an  ANOVA  was  used  to 
examine  the  variance  of  the  mean  values  for  dependent  variables  resulting  from 
different  levels  of  each  of  the  experimental  factors.  The  model  used  for  a  three-factor 
split-plot  ANOVA  for  the  dependent  variable  (Xljkl)  is  defined  below,  as  adapted  from 
Kirk,  1982. 

XtfH-ji>«J+iciW+P4+(aP)/l+(pii)BW+Yi+(«Y)Ji<+(Yn)flW+(PY)H+(apY)A,+(PY*)4aW+€^(16) 

when:  p  ljkl=  mean  estimate  of  X 

aj>Pk,Y,=  main  effects  parameters 
;ci(j)=confounding  interaction  error  parameter 
two-factor  interactions 
(P7C)k.o)’(Y^)i,(J).~  two-factor  error  parameters 
(aPy)J|d=  three-factor  interactions 
(3Y7i)kll(J)=  three-factor  error  parameter 
8jjk,=  random  effects  error 
and:  i  =  number  of  replication  =I..4 

j  =  values  for  PLANT=A,  V,  T 
k  =  values  for  %RCF=105,  125,  150 
1  =  values  for  %ARCF=0,  25,  50 

Our  null  hypothesis  (Hq)  was  that  the  values  of  pltk  are  equal  in  the  ANOVA  model 
Our  alternate  hypothesis  (HA)  was  that  the  pijk  terms  are  different. 
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The  ANOVA  of  the  split-plot  model  accounts  for  the  dependence  introduced  by  the 
presence  of  the  confounding  variable  (replications)  and  the  resulting  error  parameters 
shown  above.  However,  an  analysis  of  differences  in  the  results  generated  by  our 
algorithm  and  those  generated  by  DISASTER™  among  the  levels  of  experimental 
factors  could  not  be  examined  through  the  split-plot  ANOVA.  This  comparison  would 
have  to  account  for  another  set  of  treatments  introduced  by  the  two  solution 
procedures,  creating  an  environment  that  must  be  addressed  by  a  split-split  plot 
ANOVA  (Reynolds,  1993).  The  execution  this  analysis  is  highly  complex,  and  is 
beyond  the  scope  of  this  research  effort.  Consequently,  ANOVA  will  be  used  when 
addressing  research  objective  four,  but  not  research  objective  three. 

Assumptions  of  the  ANOVA  Model.  There  are  three  common  assumptions  for 
standard  ANOVA  models:  independence  among  groups,  normality  of  the  sampling 
distributions  of  sample  means,  and  homogeneity  of  variance.  As  previously 
mentioned,  the  condition  of  independence  among  groups  does  not  apply  in 
experimental  situations  in  which  different  treatments  are  applied  to  the  same  subject  (a 
split-plot  design)  (Cody  and  Smith,  1991:  136).  The  condition  of  normality  is  easily 
verified  through  the  use  of  a  Wilk-Shapiro  test.  All  of  our  sampling  distributions 
demonstrated  a  Wilk-Shapiro  statistic  of  greater  than  0.9,  indicating  good  normality 
(Reynolds,  1993).  All  of  our  sample  distributions  also  exhibit  good  homogeneity  of 
variance  when  examined  via  a  scatierplot. 

Post-ANOVA  Analysis.  Statistically  significant  differences  among  the  means 
associated  with  different  levels  of  the  three  experimental  factors  were  investigated 
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using  a  Tukey  Multiple  Comparisons  Method  (p=0.05).  This  test  is  used  to  rank  the 
means  and  test  for  pair-wise  differences  among  the  means  (Freund  and  Littell, 
1981:57).  Sample  means  that  differ  by  more  than  the  critical  value  (calculated  using 
the  appropriate  degrees  of  freedom  and  levels  of  significance)  are  found  to  be 
statistically  different  from  the  others  (Freund  and  Littell,  1981:56). 

The  ANOVA  also  identifies  interaction  among  the  experimental  factors.  Any 
interaction  with  a  p-value  greater  than  0.05  was  examined  by  plotting  the  means 
associated  with  each  level  of  the  experimental  factors.  The  interaction  plots  were 
visually  inspected  to  locate  cross-over  points  where  the  values  of  one  experimental 
factor  demonstrated  different  trends  given  the  values  of  the  other  experimental  factors. 
This  plot  makes  it  possible  to  identify  the  response  of  the  dependent  variable  to 
combinations  of  the  experimental  factors  that  may  have  been  nullified  in  the  ANOVA. 

Summary 

In  this  chapter,  we  examined  our  basic  algorithm  for  scheduling  a  job  shop  to 
minimize  maximum  tardiness.  This  algorithm  required  several  modifications  to 
incorporate  TOC  and  allow  valid  comparisons  with  DISA  STER™.  We  reviewed  these 
modifications,  along  with  our  experimental  design  and  data  collection  procedures. 
Finally,  we  presented  our  method  for  statistical  analysis. 
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IV  Results 


We  used  several  statistical  comparisons  to  analyze  die  results  of  our  research 
effort.  The  primary  purpose  of  these  comparisons  was  to  address  our  third  and  fourth 
research  objectives  by  comparing  the  maximum  tardiness  of  schedules  created  by 
DISASTER™  with  those  created  by  a  simultaneous  scheduling  method,  and 
determining  the  extent  to  which  the  performance  of  a  simultaneous  scheduling  method 
is  related  to  characteristics  of  the  job  shop.  The  data  used  for  these  analyses  can  be 
found  in  tabular  form  in  Appendices  C  and  D. 

Comparison  with  DISASTER™  Results 

In  order  to  asses  the  maximum  tardiness  of  schedules  created  by  DISA  STER™ 
with  those  created  by  a  simultaneous  scheduling  method  (research  objective  three),  we 
conducted  three  sets  of  comparisons.  The  data  groups  used  for  these  comparisons  are 
outlined  in  Table  4.  Details  of  each  comparison  follow  the  table. 
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Table  4:  Data  Grouping  for  Comparison  of  Schedules 


GROUP 

DEFINITION 

ALLCASE 

Best  solution  (optimal  or  heuristic)  found  by  our  algorithm  within 
the  five  minute  time  limit  for  maximum  number  of  days  tardy  for 
each  problem 

OPTCASE 

Optimal  solution  found  by  our  algorithm  for  maximum  number  of 
days  tardy  for  each  problem 

DISASTER^, 

Best  solution  for  maximum  number  of  days  tardy  for  each  problem 
as  determined  by  DISASTER™ 

DISASTER^ 

Worst  solution  for  maximum  number  of  days  tardy  for  each 
problem  as  determined  by  DISASTER ™ 

DISASTER^ 

Best  solution  for  maximum  number  of  days  tardy  for  each  problem 
as  determined  by  DISASTER ™  (excluding  problems  without 
constraint  interaction) 

%ATARDYbMl 

Percent  reduction  in  maximum  tardiness  for  each  problem  between 
our  best  solution  (optimal  or  heuristic)  and  DISASTER™'s  best 
solution  (with  or  without  interactive  constraints) 

(DISASTER^-ALLCASE) 

%ATARDYbta  =  A - - -  (  17  ) 

DISASTER 

%ATARDYWOTt 

Percent  reduction  in  maximum  tardiness  for  each  problem  between 
our  best  solution  (optimal  or  heuristic)  and  DISASTER™'s  worst 
solution  (with  or  without  interactive  constraints) 

( DISASTERman,-ALLCASE) 

%ATARDYmt  =  A - - -  (18) 

DISASTER M 

%ATARDYint 

Percent  reduction  in  maximum  tardiness  for  each  problem 
between  our  optimal  solution  and  DISA  STER™’ s  best  solution 
with  interactive  constraints. 

(DISASTER^  -ALLCASE) 

%ATARDYlm  =  A - — - -  (19) 

m  DISASTER^ 

Comparison  with  DISASTER™  Best  Case.  This  comparison  was  made  between 
the  tardiness  for  each  problem  given  the  best  solutions  (optimal  or  heuristic)  obtained 
using  our  algorithm  within  the  five  minute  time  limit  (ALLCASE)  and  the  tardiness 
for  the  same  problems  given  the  best  solutions  obtained  for  the  same  problems  using 
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DISASTER ™  (DISASTER*,^.  This  comparison  gives  some  indication  of  the  degree 
to  which  DISA  STER™  is  capable  of  achieving  the  minimum  possible  maximum 
tardiness  for  a  given  problem.  Descriptive  statistics  for  each  dependent  variable,  as 
well  as  the  percent  difference  between  the  dependent  variables  (%ATARDYbesl)  are 
provided  in  Table  5. 


Table  5:  Descriptive  Statistics  -  ALLCASE  vs  DISASTER^, 


|  VARIABLE 

MEAN  TARDINESS 

MEDIAN  TARDINESS 

STANDARD  DEVIATION 

ALLCASE 

7.2  days 

7.0  days 

2  9  days 

DISASTER** 

7.9  days 

7.0  days 

2.8  days 

%ATARDY** 

2.2  % 

0.0  % 

10.2  % 

A  histogram  of  the  percent  differences  (%ATARDYbeit)  is  shown  in  Figure  8  Positive 


values  reflect  superior  results  by  simultaneous  scheduling  of  the  constraints.  Negative 


Figure  8:  Histogram  of  %ATARDYbest 
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values  do  not  represent  superior  performance  by  DISA  STER ™  In  these  cases, 

DISA  STER ™  identified  only  one  constraint,  and  therefore  solved  a  simpler  problem 
than  the  one  solved  by  our  algorithm. 

Inspection  of  the  histogram  revealed  that  although  a  one  case  resulted  m  a 
difference  in  minimized  maximum  tardiness  as  large  of  50%,  the  majority  of  cases 
(76%)  resulted  in  no  difference.  This  observation,  along  with  the  descriptive  statistics 
indicated  that  DISASTER^  typically  performed  as  well  as  ALLCASE. 

Comparison  with  DISASTER ™  Worst  Case  This  comparison  was  made  between 
the  tardiness  for  each  problem  given  the  best  solutions  (optimal  or  heuristic)  obtained 
using  our  algorithm  within  the  five  minute  time  limit  (ALLCASE)  and  the  tardiness 
for  the  same  problems  given  the  worst  solutions  obtained  for  the  same  problems  using 
DISASTER ™  (DISASTER^J.  This  comparison  gives  some  indication  of  the  degree 
to  which  a  user  of  DISA  STER ™  may  fail  to  achieve  the  minimum  possible  maximum 
tardiness  for  a  given  problem  if  they  fail  to  schedule  constraints  in  the  best  sequence. 
Descriptive  statistics  for  each  dependent  variable,  as  well  as  the  percent  difference 
between  the  dependent  variables  (%ATARDYwom)  are  provided  in  Table  6. 


Table  6:  Descriptive  Statistics  -  ALLCASE  vs  DISASTER^, 

VARIABLE 

MEAN  TARDINESS 

MEDIAN  TARDINESS 

STANDARD  DEVIATION  | 

ALLCASE 

7.2  davs 

.7  0  davs 

2.9  days  | 

DISASTER^, 

9.1  days 

9.0  days 

3.5  days  | 

“/oATARDY.^ 

19  3  % 

15.5  % 

18.0%  | 

A  histogram  of  the  percent  differences  (%ATARDYwom)  is  shown  in  Figure  9. 
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Figure  9:  Histogram  of  %ATARDYWW11 

Inspection  of  the  histogram  revealed  a  difference  in  the  majority  (70%)  of  cases. 
For  many  cases  this  difference  was  great,  as  high  as  60%.  This  observanon,  along 
with  the  descriptive  statistics  indicated  that  ALLCASE  performance  was  typically  and 
substantially  better  than  that  of  DISASTER,,**. 

Companson  of  Optimal  and  Valid  Cases.  Out  of  the  108  benchmark  problems 
provided  by  Captains  James  and  Mediate,  39  did  not  provide  a  satisfactory  solution  for 
comparison.  In  these  39  problems,  the  best  case  solution  provided  by  DISASTER™ 
failed  to  account  for  the  existeftce  of  the  second  constraint,  thereby  eliminating 
constraint  interaction.  Additionally,  our  algorithm  was  unable  to  fmd  and  verify  an 
optimal  solution  for  26  of  the  benchmark  problems.  Therefore,  we  limited  our  third 
comparison  to  the  SI  problems  for  which  both  an  optimal  solution  was  obtained  using 
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our  algorithm  (OPTCASE),  and  interactive  constraints  existed  in  the  solution  provided 
by  DISASTER ™  (DISASTER^.  This  comparison  identifies  the  differences  in  the 
best  solutions  provided  by  DISA  STER™  and  the  optimal  value  for  minimized 
maximum  tardiness  for  a  given  problem,  without  any  of  die  complicating  factors 
mentioned  above.  Descriptive  statistics  for  each  dependent  variable,  as  well  as  the 
percent  difference  between  the  dependent  variables  (%ATARDYmt)  are  provided  in 
Table  7.  A  histogram  of  percent  differences  (%ATARDYmt)  is  shown  in  Figure  10. 


Table  7:  Descriptive  Statistics  -  OPTCASE  vs  DISASTER^ 


VARIABLE 

MEAN  TARDINESS 

MEDIAN  TARDINESS 

STANDARD  DEVIATION 

OPTCASE 

7.2  days 

7.0  days 

2.9  days 

DISASTER^, 

7.9  days 

8.0  days 

31  days 

%ATARDY„, 

2.9  % 

0.0  % 

6.5  % 

-60  -50  -40  -30  -20  -10  0  10  20  30  40  50  60 

%  Dell*  TARDY  int 


Figure  10:  Histogram  of  %ATARDYm, 
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Inspection  of  the  histogram  revealed  that  the  majority  of  cases  (71%)  resulted  in  no 
difference.  This  observation,  along  with  the  descriptive  statistics  indicated  that 
although  DISASTER ™  did  not  always  achieve  the  optimal  solution  for  minimized 
maximum  tardiness,  there  were  only  occasional  moderate  differences  between  the 
values  for  this  performance  measure  as  represented  by  the  variables  OPTCASE  and 
DISASTER^,. 

Investigation  of  A  Igorithm  Performance 

In  order  to  determine  the  extent  to  which  the  performance  of  a  simultaneous 
scheduling  method  is  related  to  the  characteristics  of  the  job  shop  (research  objective 
four),  we  examined  the  relationship  between  the  performance  of  our  algorithm  and  the 
characteristics  of  the  job  shop  designated  as  factors  in  our  experimental  design.  These 
factors  are  listed  in  Table  8. 


Table  8:  Experimental  Factors 


FACTOR 

DEFINITION 

LEVELS  1 

PLANT 

Job  shop  configuration 

A,  T,  V  | 

%RCF 

Percent  resource  criticality  factor 

105,  115,  125  I 

%ARCF 

Difference  in  percent  resource  criticality  factor 

0,  25,  50  I 

For  each  problem,  we  recorded  two  parameters  which  characterized  the 
performance  of  our  algorithm.  First,  we  recorded  the  time  required  to  find  the  best 
solution  (without  verification)  for  our  objective  of  minimized  maximum  tardiness 
(Titnd)-  This  measure  gives  an  indication  of  the  strength  of  our  branching  technique. 
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Next,  we  recorded  the  time  required  to  verify  the  best  solution  (exclusive  of  the  time 
to  find)  for  our  objective  of  minimized  maximum  tardiness  (T^^y)  This  measure 
gives  an  indication  of  the  strength  of  our  lower  bound.  If  the  best  solution  was  not 
verified  within  five  minutes  (300  seconds),  then  verification  time  was  designated  as 

^vEMFr  ~  300  ~  Tpjf/p  (20) 

For  the  26  problems  where  the  best  solution  was  not  verified  as  optimal  within  five 
minutes,  TVERIFr  represents  an  optimistic  measure  of  the  actual  time  required  to  verify 
our  best  solution.  Therefore,  the  statistical  significance  of  any  differences  in  our 
values  for  TVERlFr  across  the  various  levels  of  die  experimental  factors  are  likely  to  be 
understated  in  this  analysis. 

We  analyzed  both  dependent  parameters  with  respect  to  the  different  levels  of  the 
factors  of  die  experimental  design  using  a  split-plot  ANOVA.  The  ANOVA  indicated 
no  in  difference  TFINn  across  the  different  levels  of  the  experimental  factors. 

However,  die  ANOVA  did  indicate  a  significant  difference  in  Ty^^y  across  different 
levels  of  PLANT  (P-value  =  0.0004).  A  Tukey  test  showed  that  our  algorithm  verified 
the  optimality  of  the  best  solution  for  our  objective  function  in  problems  with  PLANT 
=  V  significandy  faster  than  in  problems  with  PLANT  =  A  or  T.  There  were  no 
significant  differences  across  dr  Terent  levels  of  the  other  two  experimental  factors. 
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Summary 

In  this  chapter,  we  presented  the  results  of  our  analyses  in  two  parts.  First,  we 
compared  the  schedules  generated  by  our  algorithm  with  those  generated  by 
DISASTER ™.  We  demonstrated  that  with  proper  constraint  sequencing, 

DISASTER ™  is  capable  of  achieving  a  level  of  maximum  tardiness  similar  to  that 
achieved  by  our  algorithm.  However,  we  also  showed  that  if  the  constraint  sequencing 
in  DISASTER ™  is  not  chosen  carefully,  our  algorithm  was  able  to  achieve  lower 
maximum  tardiness. 

Next,  we  investigated  the  solution  characteristics  of  our  algorithm.  We  found  that 
there  was  no  difference  in  time  to  find  the  best  solution  across  the  different  levels  of 
the  experimental  factors,  but  that  there  was  a  significant  difference  in  time  to  verify 
the  best  solution  across  different  values  of  PLANT,  with  V-Plant  problems  being 
verified  significantly  faster  than  T-Plant  or  A-Plant  problems. 


Table  9:  Summary  of  Results 


COMPARISON 

RESULTS 

ALLCASE  vs  DISASTER^ 

Mean  difference  =  2.2% 

Median  difference  =  0.0% 

Standard  deviation  of  difference  =  10.2% 

I  ALLCASE  vs 

DISASTER,^ 

Mean  difference  =  19.3% 

Median  difference  =  15.5% 

Standard  deviation  of  difference  -  18.0% 

OPTCASE  vs  DISASTER, 

Mean  difference  =  2.9% 

Median  difference  =  0.0%  1 

Standard  deviation  of  difference  -  6.5%  | 

TfiND 

No  significant  difference  across  levels  of  the  experimental  | 
factors  | 

f  VERIFY 

Faster  (P-value  =  0.0004)  verification  in  problems  with 
PLANT  =  V  vs  problems  with  PLANT  =  A  or  V 
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V.  Conclusions  and  Recommendcttions 


In  the  previous  chapter,  we  presented  the  results  of  our  research.  In  this  chapter 
we  will  examine  these  findings  as  they  relate  to  our  research  objectives,  and  offer 
some  suggestions  for  further  research. 

Problem  Formulation  and  Solution  Methodology 

We  formulated  the  job  shop  problem  for  die  objective  of  minimized  maximum 
tardiness  in  a  manner  consistent  with  the  principles  of  TOC.  This  formulation  differed 
from  traditional  job  shop  formulations  primarily  due  to  the  introduction  of  rods.  The 
presence  of  rods  was  reflected  by  dynamic  ready  times  among  the  constraint 
operations.  We  were  unable  to  identify  an  optimal  solution  technique  consistent  with 
our  formulation  during  our  review  of  the  literature.  We  therefore  modified  the  Graph- 
Theoretical  method  of  Florian,  Trepant,  and  McMahon  to  reflect  the  performance 
measure  of  minimized  maximum  tardiness  and  the  relevant  concepts  of  TOC. 

Comparison  with  DISASTER ™  Results 

In  general,  our  algorithm  demonstrated  little  improvement  over  DISA  STER™ s 
best  performance  for  the  objective  of  minimizing  maximum  tardiness.  A  mean 
reduction  in  tardiness  of  2.9%  was  achieved  between  DISASTER™1  s  best  solutions  for 
problems  with  interacting  constraints  and  our  optimal  solutions  for  the  same  problems. 
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In  71%  of  these  problems,  there  was  little  or  no  difference  between  the  maximum 
tardiness  of  schedules  generated  by  our  algorithm  and  those  created  by  DISASTER ™ 
This  improvement  was  less  than  we  initially  expected  since  our  algorithm  was 
designed  solely  to  rind  an  optimal  schedule  to  minimize  maximum  tardiness,  while 
DISASTER™' s  heuristic  solution  technique  attempts  to  meet  numerous  scheduling 
objectives.  Our  results  indicate  that  DISASTER™ s  heuristic  solution  technique  is 
very  strong  for  this  performance  measure. 

Despite  the  excellent  performance  of  DISA  STER™s  best  solutions  there  is  still 
room  for  improvement.  In  29%  of  die  problems  comparing  our  optimal  solutions  and 
DISASTER™ s  solutions  with  interacting  constraints,  we  achieved  at  least  10%  better 
performance  for  the  objective  of  minimizing  maximum  tardiness.  The  improvement 
demonstrates  that  simultaneous  scheduling  of  the  constraints  can  result  in  better 
schedule  performance,  which  may  translate  into  increased  system  throughput. 

Our  algorithm  demonstrated  significant  improvement  (mean  =  19.3%)  over 
DISA  STER™s  worst  performance  even  when  heuristically  limited  to  five  minutes  of 
processing.  It  was  also  observed  that  the  sequence  of  constraints  resulting  in  the  worst 
solution  for  minimized  maximum  tardiness  was  occasionally  the  sequence  which  was 
recommended  by  DISASTER™.  Our  improvement  over  DISASTER ™'s  worst 
performance  suggests  that  users  of  DISASTER™  should  examine  the  results  of  all 
possible  constraint  sequences  to  find  the  best  schedule  that  DISASTER™  is  capable  of 
producing  for  a  given  performance  measure. 


58 


Investigation  of  A  Igorithm  Performance 

We  determined  that  there  was  no  difference  across  the  various  levels  of  the 
experimental  factors  in  the  time  required  to  find  the  best  solution  for  our  objective  of 
minimized  maximum  tardiness.  However,  we  did  find  a  statistically  significant 
difference  across  levels  of  the  experimental  factors  in  the  time  required  to  verify  the 
best  solution  for  our  performance  objective  Solutions  to  problems  in  V-plants  were 
verified  significantly  faster  than  solutions  to  problems  in  T-plants  and  A-plants 
There  are  a  greater  number  of  operations  in  A-plants  and  T-plants  than  V-plants 
In  order  to  maintain  comparable  total  constraint  loading,  the  problems  with  more 
operations  have  shorter  processing  times  per  operation.  The  rods  create  a  greater 
interval  between  operations  as  the  processing  time  per  operation  decreases.  Our  bound 
does  not  account  for  the  presence  of  rods.  Thus,  it  appears  that  our  bound  is  weaker 
in  plants  with  shorter  processing  times  because  the  best  achievable  tardiness  is 
significantly  greater  (due  to  the  gaps  in  the  schedule)  than  that  estimated  by  simple 
early-due-date  sequencing,  which  ignores  the  presence  of  rods. 

Recommendations  For  Further  Research 

Scheduling  practices  using  TOC  have  not  been  subjected  to  a  great  deal  of 
empirical  research.  As  such,  there  has  been  little  information  available  to  evaluate  the 
results  of  this  process  in  a  controlled  environment.  This  research  has  demonstrated  the 
relative  strength  of  DISASTERP^s  heuristic  constraint  scheduling  technique  for  the 
objective  of  minimized  maximum  tardiness.  Continued  evaluation  of  this  heuristic  is 


59 


encouraged,  since  this  portion  of  the  software  is  a  crucial  to  determining  the 
throughput  of  the  production  system.  Any  increase  in  throughput  as  a  result  of 
improvements  to  this  heuristic  or  other  aspects  of  DISASTER ™  will  be  very  valuable 
to  users  of  this  software. 

There  are  three  general  areas  where  we  identified  a  need  for  further  research. 

First,  there  are  several  areas  in  which  our  algorithm  could  be  improved.  A  new  lower 
bound  could  be  developed  which  is  able  to  consider  the  effects  of  rods.  The  utility  of 
the  algorithm  could  also  be  enhanced  by  addressing  other  parts  of  TOC  scheduling 
such  as  subordination  of  non-constraints  to  the  constraint  schedule.  Also,  heuristic 
solution  techniques  could  be  developed  which  avoid  lengthy  optimality  verification. 
These  improvements  would  create  a  viable  scheduling  alterative  to  DISASTER™.  This 
alternative  would  not  be  dependent  on  proper  constraint  sequencing  for  maximum 
performance. 

Another  area  rich  in  opportunity  is  further  evaluation  of  DISA  STER™. 

Algorithms  should  be  developed  to  evaluate  the  ability  of  DISASTER ™  to  maximize 
or  minimize  another  objective  -  number  of  late  jobs,  makespan,  etc.  This  would  reflect 
the  other  management  concerns  relevant  to  job  shop  scheduling.  Additionally,  it 
would  be  very  useful  to  develop  an  ability  to  predict  the  small  percentage  of  cases 
where  DISASTER™'s  solutions  are  far  from  optimal  and  devise  heuristic  adjustment 
procedures  to  deal  with  those  cases. 

The  evaluation  of  DISASTER™  could  also  be  enhanced  by  the  development  of  a 
new  experimental  design  with  compete  independence  among  the  benchmark  problems. 
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A  more  thorough  analysis  could  then  be  conducted  to  determine  the  relationship 
between  the  optimality  of  DISASTER™' s  solutions  and  the  characteristics  of  the  job 
shop. 

The  final  area  of  further  research  opportunity  which  we  envision  is  in  the 
application  of  TOC  scheduling  to  the  USAF.  Multiple  constraints  are  known  to  exist 
in  many  USAF  operations.  The  budget-cutting  activities  of  1990's  can  be  expected  to 
put  more  demand  on  our  remaining  resources,  which  will  increase  the  likelihood  of 
encountering  system  constraints.  An  empirical  study  of  USAF  systems  is  needed  to 
identify  the  number  and  nature  of  existing  constraints  in  USAF  systems  (i.e.  activities 
at  an  Air  Force  Logistics  Center).  This  research  will  help  focus  the  efforts  of  both 
future  users  and  researchers  of  TOC. 

Conclusion 

This  research  effort  has  provided  the  first  empirical  evidence  of  the  effectiveness 
of  the  DISASTER™  scheduling  program.  The  existing  partnerships  with  the  Avraham 
Y.  Goldratt  Institute  will  allow  the  Air  Force  to  become  an  influential  participant  in 
the  further  refinement  of  DISA  STER™.  We  can  expect  that  the  management  of  the 
constraints  in  USAF  operations  will  contribute  to  future  success  in  executing  our 
global  mission.  As  the  Air  Force  is  challenged  to  increase  its  effectiveness  despite 
shrinking  resources,  the  strategies  expressed  in  TOC  will  be  a  valuable  source  of 
guidance  for  continuous  improvement  efforts. 
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A ppendix  A :  Job  Shop  Scheduling  Problems 


Table  10:  A-Plant  Constraint  Loading  (James  and  Mediate,  1993:  115-123) 


%RCF 

%ARCF 

Blue  Processing  Time 

Gold  Processing  Time 

105 

0 

170 

170 

105 

25 

170 

210 

105 

50 

170 

250 

115 

0 

185 

185 

115 

25 

185 

230 

115 

50 

185 

275 

125 

0 

200 

200 

125 

25 

200 

250 

125 

50 

200 

300 

Table  II:  A-Plant  Job  Order  Due  Dates  (James  and  Mediate,  1993:  115-123) 


Rep  One 

Rep  Two 

Rep  Three 

Rep  Four 

Product 

Due  Date 

Product 

Due  Date 

Product 

Due  Date 

Product 

Due  Date 

FG-D 

04  Oct  93 

FG-D 

04  Oct  93 

FG-D 

05  Oct  93 

FG-D 

05  Oct  93 

FG-D 

05  Oct  93 

FG-D 

07  Oct  93 

FG-D 

05  Oct  93 

FG-D 

06  Oct  93 

FG-D 

05  Oct  93 

FG-D 

07  Oct  93 

FG-D 

07  Oct  93 

FG-D 

07  Oct  93 

FG-D 

06  Oct  93 

FG-D 

08  Oct  93 

FG-D 

08  Oct  93 

FG-D 

08  Oct  93 

FG-D 

07  Oct  93 

FG-D 

12  Oct  93 

FG-D 

11  Oct  93 

FG-D 

12  Oct  93 

FG-D 

11  Oct  93 

FG-D 

13  Oct  93 

FG-D 

12  Oct  93 

FG-D 

13  Oct  93 

FG-D 

13  Oct  93 

mm 

14  Oct  93 

FG-D 

13  Oct  93 

FG-D 

13  Oct  93 

FG-D 

13  Oct  93 

FG-D 

14  Oct  93 

FG-D 

13  Oct  93 

FG-D 

15  Oct  93 

FG-D 

14  Oct  93 

FG-D 

15  Oct  93 

FG-D 

14  Oct  93 

FG-D 

15  Oct  93 

FG-D 

15  Oct  93 

FG-D 

15  Oct  93 

FG-D 

15  Oct  93 

FG-D 

15  Oct  93 
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Table  12:  T-Plant  Constraint  Loading  (James  and  Mediate,  1993  124-132) 


%RCF 

%ARCF 

Blue  Processing  Time 

Gold  Processing  Time 

105 

0 

360 

360 

105 

25 

360 

450 

105 

50 

360 

540 

115 

0 

400 

400 

115 

25 

400 

500 

115 

50 

400 

600 

125 

0 

430 

430 

125 

25 

430 

540 

125 

50 

430 

640 

Table  13:  T-Plant  Job  Order  Due  Dates  (James  and  Mediate,  1993:  124-132) 


Rep 

One 

Rep 

Two 

Rep 

Three 

Rep 

Four 

Product 

Due  Date 

Product 

Due  Date 

Product 

Due  Date 

Product 

Due  Date 

FG-B 

04  Oct  93 

FG-D 

04  Oct  93 

FG-D 

06  Oct  93 

FG-C 

04  Oct  93 

FG-C 

05  Oct  93 

FG-G 

05  Oct  93 

FG-F 

06  Oct  93 

FG-D 

04  Oct  93 

FG-B 

07  Oct  93 

FG-C 

05  Oct  93 

FG-D 

07  Oct  93 

FG-F 

06  Oct  93 

FG-F 

07  Oct  93 

FG-G 

06  Oct  93 

FG-F 

07  Oct  93 

FG-D 

08  Oct  93 

FG-F 

08  Oct  93 

FG-B 

07  Oct  93 

FG-C 

08  Oct  93 

FG-F 

11  Oct  93 

FG-D 

11  Oct  93 

FG-F 

07  Oct  93 

FG-B 

13  Oct  93 

FG-G 

11  Oct  93 

FG-C 

12  Oct  93 

FG-F 

08  Oct  93 

FG-G 

13  Oct  93 

FG-G 

13  Oct  93 

FG-D 

12  Oct  93 

FG-B 

11  Oct  93 

FG-B 

14  Oct  93 

FG-B 

14  Oct  93 

FG-G 

14  Oct  93 

FG-D 

11  Oct  93 

FG-G 

14  Oct  93 

FG-B 

14  Oct  93 

FG-G 

15  Oct  93 

FG-C 

15  Oct  93 

FG-C 

15  Oct  93 

FG-C 

15  Oct  93 
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Figure  13:  V-Plant  Configuration  (James  and  Mediate,  1993:133-141) 
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Table  14:  V-Plant  Constraint  Loading  (James  and  Mediate,  1993:133-141) 


|  %RCF 

%ARCF 

Blue  Processing  Time 

Gold  Processing  Tune 

105 

0 

500 

500 

105 

25 

500 

650 

105 

50 

500 

750  I 

115 

0 

550 

550  I 

115 

25 

550 

690 

115 

50 

550 

830 

125 

0 

600 

600 

125 

25 

600 

750 

125 

50 

600 

750 

Table  15:  V-Plant  Job  Order  Due  Dates  (James  and  Mediate,  1993:133-141) 


Rep  One 

Rep  Two 

Rep  Three 

Rep  Four 

Product 

Due  Date 

Product 

Due  Date 

Product 

Due  Date 

Product 

Due  Date 

FG-G 

04  Oct  93 

FG-E 

04  Oct  93 

FG-G 

04  Oct  93 

FG-A 

04  Oct  93 

FG-A 

07  Oct  93 

FG-A 

06  Oct  93 

FG-B 

05  Oct  93 

FG-A 

05  Oct  93 

FG-A 

08  Oct  93 

FG-B 

06  Oct  93 

FG-D 

06  Oct  93 

FG-G 

05  Oct  93 

FG-B 

08  Oct  93 

FG-D 

07  Oct  93 

FG-A 

12  Oct  93 

FG-B 

06  Oct  93 

FG-E 

08  Oct  93 

FG-G 

08  Oct  93 

FG-B 

12  Oct  93 

FG-E 

06  Oct  93 

FG-E 

11  Oct  93 

FG-B 

12  Oct  93 

FG-D 

12  Oct  93 

FG-G 

06  Oct  93 

FG-G 

11  Oct  93 

FG-D 

12  Oct  93 

FG-G 

12  Oct  93 

FG-B 

07  Oct  93 

FG-B 

12  Oct  93 

FG-G 

12  Oct  93 

FG-A 

15  Oct  93 

FG-E 

07  Oct  93 

FG-D 

14  Oct  93 

FG-A 

13  Oct  93 

FG-E 

15  Oct  93 

FG-D 

08  Oct  93 

FG-D 

15  Oct  93 

FG-E 

15  Oct  93 

FG-E 

15  Oct  93 

FG-D 

15  Oct  93 
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Appendix  B  Software  Code  Listing 


program  Thesis, 
uses 

WinCrt,  Win  Dos,  WinProcs; 

const 
m  =  40; 
numjob  =  10; 
bignum  =  32000; 

type 

rvector  =  array[0..61]  of  real; 
vector  =  array  [0  .61]  of  integer; 
dblvector  =  arrayf0..61,1..2]  of  integer, 
check  =  array [0..61]  of  boolean; 
matrix  ==  array[0..61,0..m]  of  integer; 
nnatrix  =  array[0..61,0..m]  of  real; 
dual  =  array[1..2]  of  real; 

var 

space,  checktime,  solvetime,  starttime,  finishtime,  mtard,  upbound,  temp: real; 
m  low  bound,  jlowbound,  lowbound,  tempbound,  nodecnt;  real; 
a,  j,  i,  n,  k,  filend,  branch,  level,  select,  ref,  dref,  genset,  rep,  termnode;  integer: 
overdue,  start,  finish:  rvector, 

job,  machine,  gensize,  duration,  due,  gold,  blue:  vector; 
pred:  dblvector; 
scheduled,  cut:  check; 

GC:  matrix; 

GCbound:  nnatrix; 

indata,  outdata:  text; 

done,  solution,  test:  boolean; 

hour,  minute,  second,  seclOO:  word; 

plant:  char, 

machfm:  dual; 


{*  NEWSTART  PROCEDURE  —Calculates  new  start  times  when  cm  operation  is  scheduled*} 
procedure  Newstart(var  nstart,  nfinish:  rvector;  var  nscheduled:  check;  var  nselect:  integer); 

var 

ment,  j,  o therm  ach  :  integer; 
tempstart:  real; 

begin 
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machfm[l]  :=  0;  { Initialize  machine  one  finish  time } 

machfin[2]  :=  0;  { Initialize  machine  two  finish  time} 

if  machine[nselect]  =  1  {Determine  machine  for  selected  operation} 

then  o therm ach  2 

else  o  therm  ach  :=  1; 

for  i  :=  1  to  n  do  { Set  finish  time  of  selected  machine } 

if  (machine[i]  =  machine[nselect])  {to  the  finish  time  of  Iciest  scheduled  operation} 
and  (nscheduledfi]  =  TRUE) 
and  (nfinish[i]  >  machfin[machine[nselect]]) 
then  machfinfmachinefnselect]]  :=  nfinish[i]; 

for  i  :=  1  to  n  do  {Set  finish  time  of  non-selected  machine} 

if  (machine[i]  =  o  therm  ach)  {to  finish  time  of  the  latest  scheduled  operation} 

and  (nscheduled[i]  =  TRUE) 
and  (nflnish[i]  >  machfmfo  therm  ach]) 
then  machfm[othermach]  :=  nfinishfi]; 

for  j  :=  1  to  2  do  {Calculate  inter-operation  spacers} 

if  pred[nselectj]  >  0  then 
begin 

space  :=  0; 

if  machine[nselect|  =  machine[pred[nselectj]| 
then 
begin 

if  duration[pred[nselectJ]]  <  duration [nselect] 
then  space  :*  (240-(0.99*duration[pred[nselectj]])) 
else  space  :=  (240-(0  99*duration[nselect])); 
if  space  <  0  then  space  :=  0; 
if  (nfmish(pred[nselectj]]+space)  >  nstart[nselect] 
then  nstart[nselect]:=(nfinish[pred[nselectj]]+space); 
end 
else 
begin 

if  duration[pred[nselect  j]]  <  duration  [nselect] 
then  space  :=  (240-(0.99*durationfpred[nselectj]])) 
else  space  :=  (240-{0.99*duration[nselect])), 
if  (nfinish[pred[nselectj]]+space)  >  nstart [nselect] 
then  nstart[nseiect]:=(nflnish[pred[nselectj]]+space); 
end; 
end; 

if  nstart[nselect]  <  machfm[machine[nselect]]  early  start/flnish  for  each  operation } 

then  nstart[nselect]  :=  machfin[machine[nselect]]; 
nfinish[nselect]  :=  nstart[nselect]  +  duration  [nselect]; 
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machfinfmachinefnselect]]  :=  nfinishfnselect]; 

nschcduled[nselect]  :=  TRUE; 

end; 


f*  MBOUND  FUNCTION  —Calculates  lower  bound  for  selected  operation  * }  function 
MBound(bselect:  integer;  bstart,  bfinish:  rvector;  bdue:  vector;  bscheduled:  check):  real; 

var 

mcnt,  complete,  short:  integer; 
done:  boolean; 

mfinish,  tardy,  tardy  max:  real; 
begin 

newstartfbstart,  bfinish,  bscheduled,  bselect);  { Calculate  new  start  times  given  selected 

operation} 

tardy  max  :=  0; 

for  i  :=  1  to  n  do  { Determine  maximum  tardiness  for  scheduled  ops} 

if  (bscheduled[i]  =  TRUE) 
and  ((bfinish[i]-bdue[i])  >  tardymax) 
then  tardymax  :=  (bfinish[i]  -  bduefi]); 

for  mcnt  :=  1  to  2  do  {For  each  machine  do...} 

begin 

mfinish  :=  machfin[mcnt]; 
done  :=  FALSE; 


while  done  =  FALSE  do 
do...} 
begin 

complete  :=  bignum; 
for  i  :=  1  to  n  do 
if  (machine[i]  =  mcnt) 
and  (bscheduled[i]  =  FALSE) 
and  (bduefi]  <  complete) 
then 
begin 

complete  :=  bdue[i]; 
short  :=  i; 
end; 

if  complete  =  bignum 
then  done  :=  TRUE 
else 


{While  due  date  calculations  are  not  done 


{Find  operation  with  earliest  due  date} 


begin 

mfinish  :=  mfinish  +  duration  [short]; 
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tardy  m finish  -  bdue[short];  {Determine  tardiness  of  selected  operation} 

if  tardy  >  tardy  max 

then  tardy  max  :=  tardy;  { Update  maximum  tardiness} 

end; 

bscheduled[short]  :=  TRUE;  {Eliminate  early  unscheduled  due  date  from 

search} 
end; 

end; 

m  bound  :=  tardy  max;  {Return  largest  machine  based  lower  bound} 

end; 


{***  GENERATE  PROCEDURE  —  Calculates  the  generating  set  for  the  current  cut  ***} 
procedure  Generate(var  glevel:  integer;  var  gcut:  check;  var  gGC:  matrix); 

var 

genmach:  integer; 
minfln:  real; 

begin 

minfin  :=  bignum;  {Select  machine  for  generating  set} 

glevel  :=  glevel  +  1;  {Increment  level  number} 

for  i  :=  1  to  n  do 
if  (gcut[i]  =  TRUE) 
and  (finish[i]<minfm) 
then 
begin 

minfin  :=  finish[i]; 
genmach  :=  machinefij; 

end; 

genset  :=  0; 
for  i  :=  1  to  n  do 
if  (gcutfi]  =  TRUE) 
and  (machine[i]  -  genmach) 
then 
begin 

genset  :=  genset  +  1; 
gGC [glevel^enset]  :=  i; 
end; 

gensize[glevel]  :=  genset;  {Save  size  of  generating  set  at  this  level} 


{Select  operations  in  cut/on  machine} 


{Increment  operations  in  generating  set} 
{Add  operation  to  generating  set} 


{Find  operation  with  earliest  finish  time} 

{ Select  machine  corresponding  to  operation} 
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if  gen  set  =  0 
then  solution  :=  TRUE 
end; 


{***  SOLVE  PROCEDURE  —  Calculates  schedule  solution  during  backtracking  mode  ***} 
procedure  solve(slevel:  integer;  sstart,  s finish:  rvector;  sblue,  sgold:  vector; 

ssolution:  boolean;  scut,  sscheduled:  check;  sGC:  matrix;  sGCbound:  nnatrix); 


var 

j,  sselect,  initbranch,  z:  integer; 
label  1; 


begin; 

solution  :=  FALSE; 
sselect  :=  sGC[slevel,branch]; 
initbranch  :=  branch; 
for  i  :=  (slevel)  to  n  do 
begin 

sblue[i]  :=  0; 
sgold[i]  :=  0; 
for  j  :=  1  to  gensize[i]  do 
begin 

sstart[sGC[ij]]  :=  0; 

sfinish[sGC[ij]]  :=  (sstart[sGC[ij]]+duration[sGC[ij]]); 
scut[sGC[ij]]  :=  FALSE; 
end; 
end; 


{ Assume  no  solution  exits} 

{Select  the  given  operation} 

{For  all  levels  below  that  of  selected  operation  do  } 
{Reinitialize  schedule  information} 


newstart(sstart,sfinish,sscheduled,sselect), fCa/cw/are  start  times} 

if  machine[sselect]  =  1  { Schedule  selected  operation  on  proper  machine} 

then  sblue[slevel]  :=  sselect 
else  sgoldjslevel]  :=  sselect; 

for  i  :=  i  to  n  do  { Select  successor  operation  for  next  cut} 

if  (sscheduled[pred[i,  1  ]]  -  TRUE) 
and  (sscheduied[pred[i,2]]  =  TRUE) 
and  (sscheduled[i]  =  FALSE) 
then  scut[i]  :=  TRUE; 

Generate(slevel,scut,sGC);  {Generate  next  Generating  Set } 

While  (solution  =  FALSE)  do  { Until  a  solution  is  reached  do...} 

begin 
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temp  :=  bignum; 


for  a:=  1  to  genset  do 

begin  {Find  the  operation  with  the  smallest  lower  bound} 

lowbound  mbound(sGC[slevel,a],sstart,sfmish,due,sscheduled); 
sGCbound[slevel,a]  :=  lowbound; 
if  lowbound  <  temp 
then 
begin 

branch  :=  a, 
temp  lowbound; 
end; 
end; 

While  temp  >=  upbound  do  {If  lower  bound  exceeds  upper  bound  then  do  } 

begin 

genset  :=  gensize[slevei];  {Reset  schedule  information  at  this  level} 

for  a  :=  1  to  genset  do 

begin 

sstart[sGC[slevel,a]]  :*  0; 

sflnish[sGC[slevel,a|]  :=  (sstart[sGC[slevel,a]]+duration[sGC[slevel,al]); 
sscheduled[sGC[slevel,a|]  :=  FALSE; 
scut[sGC[slevel,al]  :=  FALSE; 
end; 

sblue[slevel]  :=  0; 
sgold[slevel]  :»  0; 

slevel  :=  slevel  -  1;  {Backtrack  one  level} 

if  slevel  <=  level  {If  backtracked  to  subroutine  entry  point,  then  exit} 

then 
begin 

GCbound[level,initbranch]  :=  bignum; 
level  :=  level  +  1; 
goto  1; 
end 

else  { Else  try  next  lower  bound  at  previous  level} 

begin 

temp  :=  bignum; 
genset  :=  gensize[slevei]; 
for  a  :=  1  to  genset  do 
begin 

sstart[sGC[slevel,a]]  :=  0; 

sfinish[sGC[slevel,a]]  :=  (sstart[sGC[slevel,a]]+duration[sGC[slevel,a]]); 
sscheduled[sGC[slevel,a|]  :=  FALSE; 
scut[sGC[slevel,a]]  :=  FALSE; 
if  sGCbound[slevel,a]  <  temp 
then 
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begin 

branch  :=  a; 

temp  :=  sGCbound[slevel,a]; 
end; 
end; 
end; 
end; 

sGCbound[slevel, branch]  :=  bignum;  { Set  lower  bound  of  selected  operation  =  bignum} 
sselect  :=  sGC[slevel, branch];  {Selected  operation  is  to  be  schedueled  next } 

nodecnt  :=  nodecnt  +  1 ; 

if  machinefsseiect]  =  1  {Schedule  selected  operation  on  proper  machine} 

then  sblue[slevel]  :=  sselect 
else  sgold[slevel]  :=  sselect; 

newstart(sstart,sfinish,sscheduled,:sselect);/Ca/c2//dte  new  start  times} 
scut[sselect]  :=  FALSE;  {Remove  selected  operation  from  cut} 

for  i  :=  1  to  n  do  {Select  successor  operation  for  next  cut} 

if(sscheduled[pred[i,l]]  =  TRUE) 
and  (sscheduled[pred[i,2]]  =  TRUE) 
and  (sscheduled[i]  =  FALSE) 
then  scut[i]  :=  TRUE; 

Generate(slevel,scut,sGC); 

end; 

upbound  :=  0; 
for  i  :=  1  to  n  do 
begin 

start[i]  :=  sstart[i]; 
finishfi]  :=  sftnish[i]; 
blue[i]  :=  sblue[i]; 
gold[i]  :=  sgold[i]; 
cut[i]  :=  scut[i]; 
scheduled[i]  :=  sscheduled[i]; 
for  j  :=  1  to  m  do 
begin 

GC[ij]  :=  sGCfij]; 

GCbound(ij]  :=  sGCboundfij]; 
end; 

if  (sfinish[i]  -  due[i])  >  upbound  {Calculate  new  upper  bound} 
then  upbound  :=  (sfinishfi]  -  due[i]); 
end; 

GetTime(Hour,  Minute,  Second,  Secl00);fGe/  solution  time} 
solvetime  :=  ((((hour*60)+minute)*60)+second+(sec  100*0.01)); 


{Generate  new  generating  set} 


{ Update  schedule  info  with  new  schedule} 
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term  node  :=  term  node  +  1 ; 
level  :=  n+1; 

1:  end; 


{***  REINITIALIZE  PROCEDURE  -  Resets  all  variables  to  inital  status  ***} 
procedure  reinitialize; 

begin 

done  :=  FALSE; 
solution  :=  FALSE; 
level  :=  0; 
mtard  :=  0; 
nodecnt  :=  0; 
term  node  :=  0; 
upbound  :=  0; 
temp  :=  0; 
lowbound  :=  0; 
tempbound  :=  0; 
branch  :=  0; 
select  :=  0; 
genset  :=  0; 
for  i  :=  0  to  60  do 
begin 

cut[i]  :=  FALSE; 
scheduled[i]  :=  FALSE; 
job[i]  :=  0; 
machine[i]  :=  0; 
duration[i|  :=  0; 
due[i]  :=  0; 
finish [i]  :=  0; 
startfi]  :=  0; 
overdue[i]  :=  0; 
gensize[i]  :=  0; 
gold[i]  :=  0; 
blue[i]  :=  0; 
for  k  :=  1  to  2  do 
pred[ijc]  :=  0; 
for  k  :=  0  to  m  do 
begin 

GC[iJc]  :=  0; 

GCbound[iJc]  :=  0; 
end; 
end; 
end; 
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I***********************  MAIN  ROUTINE  ********************************} 
begin 

asign  (indata,'indata.txt');  {open  ascii  input  datafile} 

reset  (indata); 

assign  (outdata,  'outdata.txt');  { open  ascii  output  datafile} 

rewrite  (outdata); 

readln(indataXilend); 

while  filend  <  1  do  {While  not  end-of-file  do  .  .  .} 

begin 

gettime(hour,  minute,  second,  sec  100);  { Get  computation  start  time} 
starttime  :=  ((((hour* 60)+minute)*60)+second-Ksec  100*0. 01)); 


{***  INPUT  DATA  ***} 
read(indata,  plant,  ref,  dref,  rep,  n); 
wnteln(outdata,,================= 

writeln(outdata,plant,'-PLANT,  '/cf^/o 

wnteln(outdata); 

for  i  :=  1  to  n  do 

begin 

read  (indatajob[i]); 
read  (indata^nachine[i); 
read  (indata,pred[i,l  ],  pred[i,2]); 
read  (indata,duration[i]); 
read  (indata,due[i]); 
end; 


{Read  schedule  identification  information} 

RCF,  ',drcf,*%  DRCF,  REPLICATION:\rep); 
{  Write  schedule  header} 


{Read  the  job  number  for  each  operation} 

{ Read  the  machine  number  for  each  operation } 
{Read  the  predecessors  of  each  operation} 
{Read  the  duration  of  each  operation} 

{ Read  the  due  date  of  each  operation} 


{***  SELECT  INITIAL  CUT  ***} 
for  i  :=  1  to  n  do 
if  (pred[i,l]  =  0) 
and  (pred(i,2J  =  0) 

then  cut[i]  :=  TRUE;  { All  operations  without  predecessors  in  initial  cut} 


{***  ASSIGN  EARLIEST  START  AND  FINISH  TIMES  ***} 


start[0J  0; 
finish  [0]  :=  0; 
scheduled[0]  :=  TRUE; 


{Null  operation  starts  at  0} 
{Null  operation  finishes  at  0} 
{ Null  operation  is  scheduled} 


for  i  :=  1  to  n  do 

finish  [i]  :=  start  [i]  +  duration  [i];  { For  each  operation  finish  =  start  +  duration} 


{***  FIND  UPPER  BOUND  ***} 

generate(level,cut,GC);  {Generate  initial  generating  set} 

while  solution  =  FALSE  do  {Until  a  solution  is  found  do  for  each  level .  .  .} 

begin 

temp  :=  bignum; 
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for  a:=  1  to  genset  do 

begin  {Find  the  operation  with  the  smallest  lower  bound } 

lowbound  :=  m  bound(GC[level, a], start/mish, due, scheduled); 

GCbound[level,a]  :=  lowbound; 
if  lowbound  <  temp 
then 
begin 

branch  :=  a; 
temp  :=  lowbound; 
end; 
end; 

GCbound[Ievel,branch]  ;=  bignum;  {Remove  selected  operation  from  consideration} 
select  :=  GC[level,branch];  {Designate  selected  operation} 

nodecnt  :=  nodecnt  +  1;  {Increment  node  count} 

if  machine[select]  =  1 
then  blue[level]  :=  select 

else  gold  [level]  :=  select;  {Record  selected  operation} 

newstart(start,finish,scheduled,select);/'Ca/cu/ate  new  start  times} 
cut[select]  :=  FALSE;  {Remove  selected  operation  from  cut} 

for  i  :=  1  to  n  do 
if  (scheduled[pred[i,l][  =  TRUE) 
and  (scheduled  [pred[i,2]]  =  TRUE) 
and  (scheduled[i]  =  FALSE) 
then  cut[i]  :=  TRUE;  {Determine  next  cut} 

Generate(level,cut,GC);  {Generate  next  generating  set} 

end; 

for  i  :=  1  to  n  do 
if  (finish[i]  -  due[i])  >  upbound 

then  upbound  :=  (finishfi]  -  due[i]);  { When  solution  achieved,  calculate  upper  bound} 

gettime(hour,  minute,  second,  sec  100);  {Get  solution  time} 

solvetime  :=  ((((hour*60)+minute)*60)+second-Ksec  100*0.01)); 

term  node  :=  term  node  +  1;  {Increment  number  of  terminal  nodes} 

{***  BACKTRACKING  ROUTINE  ***} 

level  :=  level  -  1;  {Decrement  level} 

While  level  >  0  do  {While  level  greater  than  zero  do  ..  .} 

begin; 

genset  :=  gensizeflevel];  {Determine  size  of  generating  set  at  this  level} 

temp  bound  :=  upbound; 
for  i  :=  1  to  genset  do 
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scheduled[GC [level,i ] ]  :=  FALSE; 
for  i  :=  1  to  genset  do 
if  (GCbound[level,i]  <  tempbound) 
then 
begin 

branch  :=  i; 
tempbound  GCbound[level,i]; 
end; 

if  tempbound  <  upbound  then  {If  lower  bound  <upper  bound  then  goto  SOLVE} 

solve(level,  start,  finish,  blue,  gold,  solution,  cut,  scheduled,  GC,  GCbound); 
level  :=  level  -  1;  {Decrement  level} 

end; 

{***  OUTPUT  ROUTINE  ***} 

for  i  :=  1  to  n  do  {For  each  operation  do  .  .} 

begin 

overdue[i]  :=  finish [i]  -  due[i];  {Calculate  tardiness} 

if  overdue[i]  <  -240  {If  >  '6  shipping  buffer  early} 

then  overdue[i]  :=  0  {Operation  not  tardy} 

else  overduefi]  :=  overdue[i]+480;  {Operation  tardy  (add  buffer)} 
end; 

mtard  :=  0; 
for  i  :=  1  to  n  do 

then  mtard  :=  overduefi]; 

mtard  :=  mtard/480;  {Translate  maximum  tardiness  to  days} 

if  mtard  -  trunc(mtard)  >  0 

then  mtard  :=  trunc(mtard)  +  1;  {Round  up  to  nearest  whole  day} 

writelnfoutdata,'  Blue  Machine  Schedule'); 

writeIn(outdata,'operation  job  start  finish  due  tardy'); 

for  i  :=  1  to  n  do 
if  blue[i]  >  0 
then 
begin 

write(outdata,blue[i]:5);  { Output  operation  number} 

write(outdataJob[blue[i]]:9);  { Output  associated  job  number} 

write(outdata,bvart[blue[i]]:9: 1 );  {Output  operation  start  time} 

writefoutdata^nishfblueli]]^:  1);  {Output  operation  finish  time} 

write(outdata,due[blue[i]]:8);  { Output  operation  due  date} 

writeln(outdata,overdue[blue[i]]:8:l);  { Output  operation  tardiness} 

end; 

writeln(outdata); 

writelnfoutdata,’  Gold  Machine  Schedule'); 

writeln(outdata,'operation  job  start  finish  due  tardy'); 

for  i  :=  1  to  n  do 
if  gold[i]  >  0 


{Unschedule  operations  at  this  level} 

{Determ  ine  lowest  lower  bound  at  this  level} 
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then 

begin 

wnte(outdata£old[i] : 5 );  {Output  operation  number} 

write(outdataJob[gold[i]]:9);  {Output  associcted  job  number} 

write(outdata,start[gold[i]]:9: 1 );  {Output  operation  stai  time} 

wnte(outdata,fmish[gold[iJ].  8  1),  {Output  operation  finish  time} 

write(outdata,due[gold[i]]:8);  {Output  operation  due  date} 

wiiteln(outdata,overdue[gold[i]]: 8:1);  { Output  operation  tardiness} 

end; 

writeln(outdata); 

write!n(outdata,'MAXIMUM  TARDINESS  =  ',mtard:3:3,'  days') ; 
writeln(outdata); 

GetTime(Hour,  Minute,  Second,  SeclOO);  {Get  finish  time} 
fuiishtime  :=  ((((hour*60)+minute)*60)+second+(secl  00*0.0  l))-solvetime; 
solvetime  :=  solvetime-starttime;  {Calculate  solution  time} 

checktime  :=  finishtime-solvetime;  { Calculate  verification  time} 

writeln(outdata,'SOLUTION  TIME  =',soIvetime: 7:1, ’seconds'); 
writeln(outdata,'VERIFICATION  TIME  =  ',checktime:7:l,'  seconds'); 
writeln(outdata); 

write!n(outdata,'  NODES:\nodecnt:7:  l,’  TERMINAL  NODES:’,termnode); 

wnteln(outdata,'===~=================:==========='); 

writeln(outdata); 

reinitialize;  {Reinitialize  all  variables} 

readln(indata,  filend);  {Read  filend  status  variable} 

end; 
end. 
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A  ppendix  C:  Comparison  of  Solutions 
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A  ppendix  D:  TOC  Job  Shop  A  Igorithm  Statistics 
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